Quantum kinetic approach to the calculation of the Nernst effect 
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We show that the strong Nernst effect observed recently in amorphous superconducting films far 
above the critical temperature is caused by the fluctuations of the superconducting order parameter. 
We employ the quantum kinetic approach^ for the derivation of the Nernst coefficient. We present 
here the main steps of the calculation and discuss some subtle issues that we encountered while 
calculating the Nernst coefficient. In particular, we demonstrate that in the limit T — » the 
contribution of the magnetization ensures the vanishing of the Nernst signal in accordance with the 
third law of thermodynamics. We obtained a striking agreement between our theoretical calculations 
and the experimental data in a broad region of temperatures and magnetic fields. 



I. INTRODUCTION 

After many years in the shade, the Nernst effect (the 
transverse thermoelectric signal) entered the spotlight 
in condensed matter physics as well as other fields of 
research, such as the theory of gravitation^ The "re- 
discovery" of the Nernst effect by the condensed matter 
community occurred after the measurement of the effect 
in high-T c materials above the superconducting transi- 
tion temperatures^^ Since then, the Nernst effect was 
also observed in conventional amorphous superconduct- 
ing films far above T c £& The Nernst effect in high-T c 
superconductors^ has been attributed to the motion 
of vortices^^ existing even above T c (the vortex-liquid 
regime). In conventional amorphous superconducting 
films the strong Nernst signal observed deep in the nor- 
mal stated cannot be explained by the vortex-like fluc- 
tuations. Rather, the authors of Rcfs. HUD suggested that 
the effect is caused by fluctuations of the superconducting 
order parameter. Here we present a comprehensive anal- 
ysis of this mechanism using the quantum kinetic tech- 
nique and demonstrate a quantitative agreement between 
the theoretical expressions and the experiment £ No fit- 
ting parameters have been used; the values of T c and the 
diffusion coefficient were taken from independent mea- 
surements (see Refs.[l@). In particular, we succeeded in 
reproducing the non-trivial dependence of the signal on 
the magnetic field. Our results imply that in the quest for 
understanding the thermoelectric phenomena in high-T c 
materials the fluctuations of the order parameter should 
not be ignored. 

The Nernst effect and its counterpart, the Etting- 
shausen effect, are effective tools for studding the su- 
perconducting fluctuations because in metallic conduc- 
tors the contribution of the quasi-particlc excitations is 
negligible. Under the approximation of a constant den- 
sity of states at the Fermi energy, which is a standard 
approximation for the Fermi liquid theory, this contribu- 
tion vanishes completely^ On the other hand, the col- 
lective modes describing all kinds of fluctuations can in 
general generate significant contributions to the Nernst 
effect. Since the neutral modes arc not deflected by the 
Lorentz force, they do not contribute to the transverse 



thermoelectric current. The charged modes, such as fluc- 
tuations of superconducting order parameter, are a pos- 
sible source for the giant Nernst effect even far from the 
superconducting transition. The fact that the main con- 
tribution to the Nernst signal originates from the super- 
conducting fluctuations is in contrast to other transport 
phenomena such as the electric conductivity. The con- 
tributions to the electric conductivity caused by the su- 
perconducting fluctuations (paraconductivit y 11 ' 1 ^ 13 ) can 
be observed close enough to the superconducting tran- 
sition where the paraconductivity increases rapidly and 
may even overcome the Drude conductivity. Far from 
the transition the superconducting fluctuations produce 
only one among many corrections to the conductivity 
and, therefore, can hardly be identified. Owing to the 
fact that in the absence of fluctuations the Nernst effect 
is negligible, measurements of the Nernst signal provide 
a unique opportunity to study the superconducting fluc- 
tuations deep inside the normal state. 

The transport coefficients for the electric and thermal 
currents are defined via the standard conductivity tensor: 
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When the thermo-magnetic phenomena are studied in 
films (or layered conductors) the magnetic field is con- 
ventionally directed perpendicularly to the conducting 
plane, see Fig. Q] Then, each element of the conductiv- 
ity tensor corresponds to a 2 x 2 matrix describing the 
conductivity components in the x — y plane (see Fig. []}. 
The different components of the conductivity tensor are 
connected through the Onsager relations. In particular, 
(Tij(H) = H) and <Sij(H) = T<x,;(— H). In an open 



circuit setup, from the condition j e 
Nernst coefficient is: 



one gets that the 
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We checked that the second term in the numerator is 
negligible in comparison to the first one (see the com- 
This observation has been exper- 



ment below Eq. I42a| 
imentally verified as it follows from Fig. 2(a) in Rcf. 5. 
Therefore, the leading order term in the expression for 
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FIG. 1: The setup of the Nernst effect measurement. The sam- 
ple is placed between two thermal baths of different tempera- 
tures. The temperature gradient is in the z-direction, the mag- 
netic field is along the z-direction and the electric field is induced 
in the y-direction. 



the Nernst coefficient is ejv ~ ct xy /a xx and our goal is to 
find the transverse Peltier coefficient, a xy . 

The electric current generated as a response to an ex- 
ternal force, such as the electric field, can be found in 
the linear regime by the Kubo formulaM which expresses 
the response in terms of a corresponding correlation func- 
tion. Extending the Kubo formalism to the calculation 
of the response to the temperature gradient is not trivial 
because this gradient is not directly connected to any me- 
chanical force. Following the scheme used in the deriva- 
tion of the Einstein relation, Luttingcr— made a connec- 
tion between the responses to the temperature gradient 
and to an auxiliary gravitational field. As a result, Lut- 
tinger succeeded in relating all transport coefficients with 
various current-current correlation functions. A main in- 
gredient in the Kubo formula is the quantum mechanical 
expression for the current operators that enter the cor- 
relation function; e.g., the electric and heat currents in 
case of the thermoelectric transport . When the electron- 
electron interactions are neglected, the expression for the 
heat current operator is 

j h (q = 0, r) = J2 tt^p - vKJtKAt) (3) 
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dp 
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where c p CT (r) (c P (T (t)) is the creation (annihilation) op- 
erator of an electron in a state with energy e p . Here y, is 
the chemical potential, Vi mp (p, p') is the potential cre- 
ated by the disorder and r is the imaginary time. With 
the help of the equations of motion and after transform- 
ing to the Matsubara frequencies, the current operator 
can be written as>i£ 



j,,(q = 0,u> n ) = 



de p 2ie n - iu> n 
~dp~ 2 



c P ,o-( e ™) c P,<x( e n - 

(4) 



When electron-electron interactions are included, is 
more complicated function (it contains terms with four 
fcrmion operators). In general, the resulting expression 



for the heat current is not just the frequency multiplied 
by the velocity as it is for free electrons. Unfortunately, 
very often the expression for the heat current of free 
electrons presented in Eq. U is used in the presence of 
electron-electrons interactions, when there is no real jus- 
tification for it. In Appendix B of Ref. [l3 we showed 
that this simplified form of the Kubo formula fails to re- 
produce the thermal conductivity of Fermi liquids. The 
incorrect result that emerges from Eq. [4] does not imply 
that the use of the Kubo formula for the thermal trans- 
port coefficients is necessarily wrong. The weak point is 
in replacing the full expression for the heat current by 
the one in Eq. [4j The problem with the full expression 
for the heat current is in its complexity. 

In addition, the Kubo formalism meets with some dif- 
ficulties when the thermoelectric currents are considered 
in the presence of a magnetic field. Obraztsov^ pointed 
out that when a magnetic field is applied, the heat cur- 
rent describing the change in the entropy must include 
a contribution from the magnetization. This is because 
the thermodynamic expression for the heat contains the 
magnetization term. Thus, additional problem of the 
Kubo formula is that the current cannot be expressed 
entirely by a correlation function. In order to determine 
the transverse thermoelectric currents one needs to com- 
bine the quantum mechanical response to the external 
field with the magnetization, which is a thermodynamic 
quantity-±&±&22 

In the derivation of the thermoelectric currents we de- 
cided, instead of applying the Kubo formula, to employ a 
different approach and to use the quantum kinetic equa- 
tion,2Ii22,>22, One main advantage of the quantum kinetic 
approach is that the problem of the magnetization cur- 
rent is solved straightforwardly. We directly obtained 
the expression for the thermoelectric current which in- 
cludes the magnetization current. In this way, the elec- 
tric current generated by the temperature gradient can 
be related to the flow of entropy. Therefore, according to 
the third law of thermodynamics the Nernst signal must 
vanish at T -> 0.24 As we will see, this argument im- 
poses a strict constraint on the magnitude of the Peltier 
coefficient in a broad range of temperatures. Note that 
the calculation of the thermoelectric transport using the 
kinetic equation allows a direct verification of the On- 
sager relations between the off-diagonal components of 
the conductivity tensor (sec Appendix [A"]l . 

The paper is organized as follows: in Sees. ITT1 and Hill 
we present the main steps in the derivation of the elec- 
tric current as a response to a temperature gradient in 
the presence of fluctuations of the superconducting order 
parameter using the quantum kinetic equation. Then, 
in Sees. IIVI and |V] we give details of the calculation that 
are specific to the transverse current. We devote Sec. IVII 
and Appendix [B] to the contribution of the magnetiza- 
tion current to the transverse thermoelectric current. We 
demonstrate that the magnetization current ensures the 
vanishing of the Peltier coefficient in the limit T — > 0. 
This makes the Nernst signal compatible with the third 



3 



law of thermodynamics. The result of the calculation of 
the Nernst effect and comparison with the Nernst sig- 
nal measured in amorphous superconducting films^ are 
presented in Sec. IVHl The content of Sec. IVIII has al- 
ready been published 2 ^ as a separate letter; we include 
it here for completeness. In Appendix [X] we demon- 
strate that the two off-diagonal coefficients of the con- 
ductivity tensor that are found independently using the 
quantum kinetic approach satisfy the Onsager relations, 
otij (B) = T5iji{— B). In view of the frequently used argu- 
ment that the particle-hole symmetry limits the magni- 
tude of the Nernst effect (see e.g. Ref. [26h we discuss this 
issue in Appendix [C] We demonstrate that the value of 
the Nernst coefficient is not constrained by the particle- 
hole symmetry. Rather, the contribution from the quasi- 
particle excitations is zero when their density of states 
is taken to be constant, which is often confused with the 
particle-hole symmetry. 



II. THE QUANTUM KINETIC EQUATION 
ABOVE T c IN THE PRESENCE OF A 
TEMPERATURE GRADIENT 

In this paper we extend the scheme developed in 
Ref. [13 to the case of electrons interacting with supercon- 
ducting fluctuations in the presence of a magnetic field. 
Here we describe the system using two fields; one is the 
quasi-particlc field -0, while the other represents the fluc- 
tuations of the superconducting order parameter A. The 
matrix functions G(r, t; r', t') and L(r, t; r', t') written in 
the Keldysh for m 21 i 22 ' 23 describe the propagation of these 
two fields, respectively. Throughout the paper, we work 
in the basis of the retarded, advanced and Keldysh prop- 
agators: 

C*(r,t,T,t)- y G A (r,t;r',t') ) ' (5) 

where a similar expression can be written for L. [Notice 
that we use the term propagators when referring to both 
these functions, while separately we name G(r, t; r', t') 
the quasi-particle Green function and L(r, t;r',t') the 
propagator of the superconducting fluctuations.] The 
derivation of the transport coefficients in the quantum 
kinetic equation is separated into two steps. First, the 
propagators are found using the quantum kinetic equa- 
tions. Then, the expression for the current in terms of 
the propagators is derived. 

We now derive the quantum kinetic equations for the 
propagators G and L in the presence of a temperature 
gradient. Inspired by Luttinger— , we introduce an auxil- 
iary gravitational field of the form 7(1") = Tq/T{t) (where 
To is the constant part of the temperature) . The purpose 
of the gravitational field is to compensates for the non- 
uniform temperature at the initial state. In other words, 
the temperature gradient and the gravitational field are 
applied in such a way that at t — — 00 the system is in 



equilibrium. Then, starting at t — —00, the gravitational 
field is adiabatically switched off. From the response to 
switching off the gravitational field, we can learn about 
the effect of the temperature gradient on the system (for 
more details see Ref. Il7[ ) . The general expression for the 
action in the presence of a gravitational field 7(1") and a 
vector potential A(r) is 



S = / drdt^r) \ 



-i-^t( r t )|L ( P t ) 
7(r) at 



(6) 



-^V--A(r)^M) 
-^™ P (r)^^(r,t)Vv(r, t) 



-~(A(r,t)il>t(r,t)^_ a (r,t) + h.c) 



|A(r,t)p 
A 



Here A is the coupling constant of the interaction (we are 
interested in the case of an s-wave coupling). The choice 
of signs is such that A > corresponds to an attractive 
interaction. The spin index a = I (— I), or equivalcntly 
t (|), indicates the spin direction up (down). In the 
above equation and throughout the paper we set H = 1. 

The Dyson equation for the Green function in the pres- 
ence of the gravitational field is: 



8 



I 

2m 



dt 

-7(r) (V imp (r) - m) 



A(r) 7 (r) V - -A(r 



(7) 



G(r,t;r',t')=6(r-r')6(t-t') 



+ 7 (r) / dt 1 dr I E(r,*;r 1 ,* 1 )7(r 1 )G(r IJ t I ;r / ,t')- 

In general, the Green function G and self-energy E con- 
tain spin indices. Since we do not consider scattering 
mechanisms that flip the spins and ignore the Zeeman 
splitting, here and in the following we do not indicate 
the spin indices whenever it is possible. Next, we intro- 
duce the following transformation: 



Y(v, t; r', = 7 " 1/2 (r)£(r, t; r', fh^V), 



(8) 



where Y can be either G or E. For the calculation of the 
response to switching off the gravitational field in the 
linear regime, we set 7(1") = 1 + r • VT/To. Then, the 
quantum kinetic equation for the Green function of the 
quasi-particlcs becomes: 



„ r ■ VT\ d 1 /„ xe . 
\l-—^)- + — V--A(r,t) 



T / dt 2m 



(9) 



■V imp (r) + fi 



G{r,t;r',t')=6(r-r')6(t-t') 



+ J dt 1 tir 1 £(r,t;r 1 ,t 1 )fi(r 1J t I ;r',0- 

The dependence of this equation on the temperature 
gradient is much simplified by the transformation to 
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(b) 



FIG. 2: (a) Illustration of Eg. llOal for the full Green function G. 
(b) The Dyson equation for Gi„t (see Eq. I10b|l . Note that Gi n t 
includes scattering by impurities only through £(G) which is a 
function of the full Green function G. The bare Green function, 
i.e., free from the interactions and the scattering by impurities, 
is denoted by Gt- 



G and S because VT/To was eliminated from all the 
terms in the equation except the derivative with re- 
spect to time. In field theories which include a non- 
trivial space-time metric g the transformation of the kind 
YJx, x') = y/— det g{x)Y (x, x')y — det g(x') (where x is 
a 4-vector) is standard. The success of this transforma- 
tion in simplifying the quantum kinetic equation is due 
to the relatively simple structure of the metric. 

We write Eq. [9] in coordinate space because all the 
propagators and self energies that enter the kinetic equa- 
tion are not translationally invariant. There are three 
sources for the lack of translation invariance. The first 
one is because the propagators and self energies depend 
on the magnetic field through the vector potential, which 
is a function of the coordinate. The second reason is due 
to the fact that we did not yet perform the averaging 
over the disorder. Finally, and most important, in the 
presence of a temperature gradient (even in the absence 
of a magnetic field) the propagators become functions of 
the center of mass coordinate. 

We choose to postpone the averaging over impuri- 
ties until the last stage of the derivation of the current. 
Therefore, the Green function of the quasi-particles con- 
tains open impurity lines as illustrated in the two coupled 
equations presented in Fig. 



G(r,t;r',t') = G int (v,t;r',t') 



(10a) 



dr.dt.Gintir, t; r 1; t 1 )^„ lp (r 1 )G(r 1 , t x ; r', t'); 



G int (r,t;r',t') = G b (r, t; r', t') + ldr 1 dt 1 dr 2 dt 2 G b (r,t;r 1 ,t 1 ) 

x E(r 1> t 1 ;r aj t a )6 <nt (r a ,t 9 ;r' > t / ). (10b) 

Here, Gj„t(r, t; r', t') is the Green function of interacting 
electrons, while G&(r, t; r', t') is free from both the in- 
teractions and the scattering by impurities. Note, that 



Gi n t{r, t; r', t') includes partially the scattering by impu- 
rities. 

Next we write the quantum kinetic equation using the 
center of mass coordinates for space and time, R = 
(r + r')/2, T = (t + t')/2 and the relative space and time 
coordinates, p — r — r', r = t—t'. Since the gravitational 
field is independent of time and we are interested in the 
steady state solution, the Green function will be taken 
to be independent of T. On the other hand, the depen- 
dence of the Green function on R remains, because the 
temperature gradient enters the equation as the product 
r • VT = (R + p/2) • VT. This dependence on R is the 
main difference between the response to a temperature 
gradient and the response to an electric field. The point 
is that in the presence of an electric field the quantum ki- 
netic equation can be formulated in such a way that the 
electric field enters only as a product with the relative 
coordinate, (r — r') • E. Therefore, after averaging over 
the disorder the electric field dependent Green function 
becomes translationally invariant. 

In order to find the expression for the VT-dependcnt 
Green function using the quantum kinetic equation, we 
separate the Green function into three parts: 



Q- — 9eq + Gl oc - eq + GvT- 



(11) 



The first part describes the propagation at equilibrium. 
The retarded and advanced components of g eq are de- 
scribed by Eq. O with VT = 0: 



e+ J-(V--A(R + p/2,t) 
2m \ c 



(12) 



-V imp (R + p/2) + n 



- J dr 1 <T^ A {p -r„e; A, imp)g^ A (r 1: e; A, imp) = 5(p). 

This is the usual Dyson equation for the Green function 
at equilibrium, in which we performed the Fourier trans- 
form of the relative time r. In the above equation we 
introduced the equilibrium self-energy, a eq . The Green 
function g eq depends on the center of mass coordinate 
through the vector potential and the potential Vi„ lp cre- 
ated by the impurities at a specific realization. Corre- 
spondingly, we use the notation g eq (p, e; A, imp) in which 
these dependencies on R are incorporated into A and 
imp. The gradient, V = | V,r + V p , in the equation for 
gf?q contains the derivatives with respect to both R and 
P- 

According to the standard rule, the Keldysh compo- 
nent of the Green function at equilibrium can be written 
in terms of the Fermi distribution function n^e) and the 
retarded and advanced Green functions: 



9e q (Pi e; A imp) = (1 - 2n F (e)) 

x [9eq(P: e 5 A imp) - gtq(p, e ! A , imp)] 



(13) 
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In the presence of a uniform and constant in time mag- 
netic field, the expressions given in Eq.rT2]can be rewrit- 
ten as a product of the phase 



(14) 



exp 



eB 

-i — • \(r — r') x (r + r') 

4c 



and the gauge invariant Green functions, g. The retarded 
and advanced components of g satisfy the equation: 



2m 



eB x p 



Vi 



imp 



(15) 



where the product of the Green function and the self- 
energy should be understood as a convolution in real 
space (see Eq. [T"2"|) . In the following the permeability 
is taken to be 1, and correspondingly we will not dis- 
tinguish between B (the magnetic flux density) and the 
magnetic field H . After averaging over the disorder, the 



gauge invariant Green functions at equilibrium become 
translationally invariant, i.e^ functions of the relative co- 
ordinate p alone (see Ref . [271 and references therein) : 



1 

2m 



dp 2 



3 2 H 2 p 2 



Ac 2 



^2r 



r R,A 



X 9, 



~"% A (p,e) = 5(p), 



(16) 



where r is the elastic mean free time of the electrons. 

As we have already discussed, when we turn from the 
equilibrium Green function to the VT-dependcnt Green 
function, an additional dependence on the center of mass 
coordinate appears. We wish to isolate this dependence 
on R from the others. Similar to <? eg (p, e; A, imp), we 
denote the dependencies of G on the center of mass co- 
ordinate caused by the impurity potential and the vector 
potential by imp and A, respectively. Then, the remain- 
ing explicit dependence on R in G(R; p, e; A, imp) arises 
due to the temperature gradient. Therefore, in the pro- 
cess of linearizing the equation in VT/Xb, we expand G 
and £ in the collision integral with respect to this explicit 
dependence on R. In other words, we may write 



J dr-til ^R + p — t 1: e; A, imp^j G_ ^R — ^—^ — -; r ± , e; A, imp^j as J dr 1 t±(Tl;p — r 1 ,e]A,imp)G(R\r 1 ,e]A,imp) 

p r 1 9^(R;r 1; e; A, imp) 



r 1 <9£(R; p — r 1; e; A, imp) „ 
dr 1 — ■ — — G_(R; r l7 e; A, imp) - j dr^R; p - r 1 , e; A, imp) - 



an 



(17) 



r 



As we shall see below, the last two terms in the expansion 
are actually proportional to VT/T . 

The equation for the local equilibrium Green function 

is: 



dr^^ 1 (p - r x , e; A, imp) Gi oc ^ eq (R; r 1 ,e; A, imp) 

dr 1 E ioc _ e? (R; p — r 1; e; A, imp) g eq (r u e; A, imp) 
R VT 



-eg eq {p,e;A,imp). 



This equation is solved by: 



(18) 



A a ■ \ R-VT dg eq (p,e;A,imp) 
Gi oc - eq (R; p, e; A, imp) = e- 



dc 



(19) 



where the corresponding self-energy should 
be taken as E; oc _ e9 (R; p, e; A, imp) = — (R • 
VT/Xb)e<9<T e9 (p, e; A, imp)/de. We see that the lo- 
cal equilibrium Green function is a straightforward 
extension of the equilibrium Green function for a non- 
uniform temperature. Since the same holds for E; oc _ e9 , 



the equation for G\ oc - eq is a closed equation determined 
by the equilibrium properties of the system. 

The Green function Gi oc - eq describes the readjustment 
of quasi-particles to the non-uniform temperature when 
the system is trying to maintain a local equilibrium. This 
response of the electrons to the temperature gradient 
tends to induce modulation of the density. Since for 
charged particles it is impossible to have a large scale 
charge modulation, the temperature gradient transforms 
into a gradient of the electro-chemical potential. There- 
fore, j e = <t(E — V/i/e) = <tE* where the effective field 
E* is the one measured in experiments. 

The role of the local-equilibrium Green function is 
most peculiar when the response to the temperature 
gradient is considered in the presence of a magnetic 
field. Under these conditions, as we show in Sees. IIIII 
and IVII G/ oc _ eg (R; p, e; A, imp) is responsible for the 
non-vanishing contribution to the electric current from 
the magnetization. 

All the remaining terms in the quantum kinetic equa- 
tion determine the last term of the Green function, Gvt: 



6 



drig~q (p - r lf e; A, imp) G V t (R; r If e; A, imp) ^f—eg eq (p, e; A, imp) 



-L(A_!f A(R + p/2) 

zm \c»p c 



2T 

dGz oc _ eg (R; p, e; A, imp) 



(20) 



an 



= / dr-jSvT (R; P - r 1; e; A, imp) <? eg (r 1; e; A, imp) 



rj dt,ioc-eq (R;p- i"i,e; A, imp) 



2 

dr^eg (p — r 1; e; A, imp) 



5e g (r l5 e; A, imp) 



SR 

p-rj dGioc-eq (R , e; A, imp) 



<9R 



r 



In the above equation, the derivatives with re- 
spect to the center of mass coordinate act only on 
the explicit dependence of G/ oc _ eg (R, p, e; A, imp) and 
£ioc-eg(R> Pi e i A, imp) on R (i.e., through the spatially 
dependent temperature). Recall that the derivatives with 
respect to the center of mass coordinate which act on 
Vi m p and A in the local equilibrium Green function was 
already included in g~ q that appears in Eq. [T8l 

Once the explicit expressions for G\ oc _ eq and E/ oc _ eg 
are inserted, the equation becomes much simpler: 



Gvt(p, e; A, imp) = g eq (e) £ V t (e) <? e9 (e 
. VT 



(21) 



'2T 



<9<? eg (e) „ , , > . , > . , ,dg eq (e) 
v eq {e)g eq (e)-g eq {e) v e? (e)- 



(9r 



The product of matrices should be understood as a convo- 
lution in real space. The velocity v eq is the renormalized 
velocity at equilibrium: 

V — V' 

v eg (r, i; r , t) = -i lim — i(r - r')a eq {r, r , e). 



(22) 



Let us point out an important difference between the 
two parts of the Green function depending on the tem- 
perature gradient, Gi oc - eq and Gvr- As was already 
mentioned, Gi oc - eq and Ei oc _ e<? are a straightforward ex- 
tension of the equilibrium Green function and self-energy 
for a non-uniform temperature. On the other hand, the 
equation for Gvt contains the VT-dependent self-energy 
which by itself is a function of Gvt- Thus, this is a self 
consistent equation, and in order to find a close expres- 
sion for Gvt, one has to determine the structure of the 
self-energy. Once the form of the self-energy is known, 
one should take into consideration in the coarse of lin- 
earization with respect to VT that all the Green func- 
tions in Svt may depend on the temperature gradient. 

To complete the derivation of the electric current as 
a response to a temperature gradient, we must also find 
the dependence of the propagator of the superconduct- 
ing fluctuations L(r, t;r',t') on VT. In the regime of 
linear response, the explicit dependence on the tempera- 
ture gradient can be eliminated from the kinetic equation 



for L by transforming to the propagator L: 

dr 1 A- 1 (r-r 1 )|(r 1 ,t;r',t / ) (23) 
<5(r-r')- / dr 1 tft i n(r,t;r I ,t I )£(r 1 ,t I ;r , ,t / ). 



Thus, the entire dependence of the propagator on the 
temperature gradient is through the self-energy term II, 
which is a function of the quasi-particle Green functions. 

Let us separate the solution of Eq. [53] into the equi- 
librium and VT-dcpcndent propagators, L = L eq + 

Lioc-eq + Tvt- The propagator at equilibrium satisfies 
the equation: 

-l 



V eq (R;p,u) = U- 1 (p)+U eq (R;p,u) 



(24) 



The entire dependence of L eq (uj) on the frequency is 
due to dressing of the bare propagator by its self-energy 
H eq (uj). In the above equation the propagator of the 
superconducting fluctuations is a function of the temper- 
ature To. Similar to Eq. [221 we mav define the "renor- 
malized velocity" of the collective mode describing the 
superconducting fluctuations at equilibrium to be: 



V eq (r, t; r', t') = -i(r - r')n e ,(r, t; r', t' 



(25) 



Note that in fact V docs not have the dimension of a 
velocity. 

The equations for the VT-dcpcndent propagators re- 
mind the first term in Eq. [21] for Gvt: 



i; oc _ e(? (R; p, lu) = -L eq (uj)Ili oc - eq (uj)L eq (uj), (26) 



and 



T V t(R;P,^) = -L eq (uo)Xl^ T {u>)L eq {u>) 



(27) 



Once again, one should understand the product as a con- 
volution of the spatial coordinate. 



III. THE ELECTRIC CURRENT AS A 
RESPONSE TO A TEMPERATURE GRADIENT 

For the calculation of the Nernst effect we have to de- 
rive the expression for the electric current as a response 
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to a temperature gradient. In the presence of a magnetic 
field, the electric current is a sum of two terms: 



of two terms: 



' icon i imaq 
Je = J e +J e ■ 



(28) 



The first one, je° n , is derived using the continuity equa- 
tion for the electric charge. The second contribution 
to the electric current originates from the magnetization 
current. Since the magnetization current is divergence- 
less it cannot be obtained using the continuity equation 
and it will be found separately. 

As follows from the action in Eq. [6l the fields ip and A 
carry electric current. Therefore, the charge continuity 
equation must include both fields: 



eJ2dt\Mr,t)\ 2 + V- jf"(r,t) = -2 ie7 (r) (29) 

a 

[At (r, t)xP l (r, i)V> T (r, t) - A(r, t)i/>\ (r, t)^\ (r, i)] , 



where if n {v,t) = -ze £ ct 7WM(r, t)V^ a (v, t) - 
Vij)l{T,i)i)„{T,t) - 2ieA|-0 o .| 2 /c]/2m. The terms in the 
RHS describe absorption and emission of quasi-particles 
by the superconducting fluctuations; the factor 2 reflects 
the fact that the Copper pairs carry charge of 2e. 

To find the expression for j™' 1 in terms of the Green 
function, we rewrite the charge density using the lesser 
component of the Green function: 



p(v,t)=-e lim ^<^(r',t')^(r,t)) (30) 
= ie lim ^2G < (r,t;r',t'). 



r — >r 



We use the notation t' — > t + to indicate that the limit 
should be taken in such a way that t is on the upper 
branch of the Keldysh contour, while t' is on the lower 
branch.— The summation over the spin projection results 
in a factor of 2. Here (A) denotes the quantum mechan- 
ical averaging with the action given in Eq. [6] Therefore, 
the Green function is fully dressed by the interactions 
and depends on the impurity potential. In addition, G is 
a function of the temperature gradient. Since we find the 
current by extracting it from the continuity equation, we 
assume that the temperature gradient has some spatial 
modulations that will be set to zero at the end of the 
procedure. 

Next, we insert the above expression into the continu- 
ity equation given in Eq . |2"91 and rewrite V -jf n as a sum 



d d 



(31) 



2 1 = e lim 

r — >r 

t'->t+ 



at7 ]£v^(r',OVvM)) 

i lim ( 7 (r')At(r',i')Vi(r,i)^ T M) 
-^(r',t')^(r',t')7(r)A(r,0; 



r — >r 



To resolve the expression for j^ ™ we need to find the 
equations of motion for the field ip. The variational 
derivative of the action in Eq. [6] with respect to the ijr 
yields the equation of motion for the field ip: 



^^('-! A W>(<v-|A W >„(M) 

+ 7(r)tW(r)W(r, t) + <7 7 (r)A(r, t)^_ a (r, t). (32) 



Under the average, the equations of motion allow us to 
rewrite the expression for X x as: 



Zj = — ie lim ( 



zl(v-^A(r)) 7 (r)(v-<5 A (r) 



-J_(v'-^A(r')) 7 (r')(V-^A(r') 

2m V c / V c 



+j{r)V lmp {r) - j(r')V imp (r') 



(T 

£ a 7 (r')At(r', i')V-<r(r', OlM'. i) 



(33) 



Now, we wish to express the electric current in 
the presence of a gravitational field in terms of the 
propagators. The expression (A(r, t)ip^ a (r', t')^^, t)) 
and its counterpart arc averaged with respect to the 
Hamiltonian that includes both the interactions and the 
gravitational field. These expressions can be written 
in terms of the self-energy of the quasi-particles. 
For example, _ (a 7 (r)A(r, t)iPl(v', t')^_ a (v, t)) 

-i 7 (r) / dr 1 dt 1 S <T (r,t;r 1 ,t 1 ) 7 (r 1 )G CT (r 1 ,i 4 ;r',i')- (Here, 
the factor i appears because in real time the evolution 
operator is of the form e~' int and because of the conven- 
tional definition of the propagators and self-energies^). 
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As a result, we obtain: 

V7(r)(V - teA(r)/c) 



X 1 = — 2e lim 

r — >r 
t'—tt 



V'7(r')(V' + ieA(r')/c) 
2m 



2 m 



G(r,i;iV) 



G(r > t;r , ,f) 



+ 7 (r) / dr 1 (ft 1 S(r,t;r 1 ,t 1 )7(r 1 )G(r 1> t t ;r / ,t') 



dr.dt.Giv, t; ri , i 1 ) 7 (r 1 )E(r 1 , t t ; r', i')7(r') 



(34) 



The factor 2 is a consequence of the sum over the spin 
index. Similarly, we can express the averages in the equa- 
tion for X 2 in terms of the self-energy of the superconduct- 
ing fluctuations, e.g., (7(1"') A^(r', t')if>i(r, t)ip^ (r, £)) = 

-i J dr 1 di 1 fl(r,£;r 1 ,£ 1 )7(i- 1 )i(r 1 ,i t ;r',i')7(r')- 

In the regime of linear response we may eliminate the 
explicit dependence of the current on 7 (r) by expressing 
the current in terms of G, L, S and II (as defined in 
Eq.[5]). Then, the sum of the two contributions to V-jg°™ 
becomes: 



~2e lim 

r — >r 
t'— >i 



V(V - zeA(r)/c) 



2m 



G(r,t;v',t') 



V(V' + ieA(r')/c) 
2m 



G(r,t;r',t') 



dr^uMr^; T l ,t l )§L(T l ,t l -,T , ,1?) 
dr I dt 1 £(r,i;r I ,i I )£(r I ,t I ;r , J t') 
drjcftiXi^r, t; r l7 t 1 )L(r 1 , ^; r', t') 

dTidULir, t; r 1 , ti)II( r i , ^; r', t') 



(35) 



Note that the current still depends on the gravitational 
field (i.e., on the temperature gradient) through the prop- 
agators and self-energies. 

In the final step of the derivation one has to resolve the 
expression for the current out of the gradient. In other 
words, to reformulate the products of propagators and 
self-energies as a derivative with respect to the center of 
mass coordinate. As discussed in the previous section, we 
can isolate the dependencies on the center of mass coor- 
dinate created by the vector potential and by the impuri- 
ties. After we average the current over the disorder and 
transform to the gauge invariant propagators and self- 
energies, these dependencies vanish. Pay attention that 
when the limit r' —> r is taken, one may rewrite Eq. l35l in 
terms of the gauge invariant quantities alone. Therefore, 
we expand the products of the propagators and the self- 
energies with respect to the deviation from R exactly in 



the same way as performed in Eq. I17I As a consequence 
of the symmetric form of the terms in Eq.[35l [SG— GE] < 
and [I1X — LIL] < , one may check that all even orders in 
the expansion vanish. In the regime of linear response it 
is enough to keep only the lowest non- vanishing order in 
the expansion. Eventually, the expression for the current 
becomes: 

jf*(r, t)=ie( dv'dt' [v(r, t; r', t')G{v', t'; r, t)l < (36) 



h.c 



ie / dv'dt' V(r,t;r',t')L(r',t';r,t) 



We use the notation [•■■] < to remind that the expression 
inside the square brackets is a product of matrices and to 
indicate that the current corresponds to the lesser com- 
ponent of the resulting matrix. The matrices v and V 
are the renormalized velocities defined in Eqs. |2~2*1 and |2~5| 
with the VT-dependent self-energies S and II replacing 
the equilibrium ones. 

The velocity of the quasi-particles v is renormalized 
by the self-energy, <5v(r, t; r', t') = —i(r — r')S(r, t; r', t'). 
We find it useful to rewrite this expression as follows: 
-[i(r'-r) + 2i(r-r')]S(r,t;r / ,t'). The idea behind this 
representation can be explained using, as an example, the 
first order expansion of the self-energy with respect to the 
superconducting fluctuations presented in Fig. [3l In this 
example, the self-energy contains a quasi-particle Green 
function propagating from r' to r and a propagator of 
the superconducting fluctuations that goes from r to r'. 
Correspondingly, the first difference of the coordinates in 
the square brackets acts on G, while the second (which 
appears with the factor 2) acts on L. 

As mentioned in the beginning of this section, there is 
another contribution to the electric current arising from 
the magnetization: 



•mag 



-2icV x M(r) lim G(r' ,t';r,t) 



(37) 



where M(r) = er x v/2mc denotes the magnetization 
and the factor of 2 is due to the summation over the spin 
index. We would like to emphasize that since the mag- 
netization is created by itinerant electrons, the magne- 
tization current can equally contribute to the transverse 
transport electric current. 



IV. DERIVATION OF THE TRANSVERSE 
COMPONENT OF f° n 

At this stage of the derivation we shall consider only 
jcon jjggpjng f or a w hile the magnetization current aside. 
Inserting the expressions for the VT-dependent propaga- 
tors given in Eqs. I2T1 [T^l and |2"T1 into Eq . I3"6l and extract- 
ing the lesser component, we obtain )f° n as a response to 
the temperature gradient. First of all, one may observe 
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that the contributions of the local equilibrium functions terms that vanish after averaging over the volume, the 
Gi oc -eq and Li oc _ eq to vanish upon averaging the non-zero part of jg° n becomes: 
current over the sample. Since we are not interested in 



if? = I K(e) 9 R (e)v?(e)g A (e) + vf{e)g%)vf{e)g%) - v R {e)g R {e)vf(e)g R {e) (38) 



27! 



2tt de 



-g R (e)v R (e)g R (e)v A (e)] J ^en F (e) 

de 



V R (e)^-v R (e)g R (e) v R {e )g R (e)v R (e)^ 



~ le j 2^ (e)ffi?(£) [ E vtW(1 - n F {e)) + Z> T (e)n F (e)] (g R (e) - /(e)) 
+ m l f^H^M [ n vT(^)(l + n P {w)) - n> T (Lu)n P (uj)] (L r {oj) - L a {oj)) + c.c. 



Here and from now on we omit the notation eq from 
the equilibrium quantities such as the propagators, self- 
energies and velocities. 

As we are interested in the Gaussian fluctuations, we 
replace the equilibrium Green function by g(r, r',e) = 
<7 (r,r',e) + / dr^go^, r 1; e)<r(r 1 , r 2 , e) g (r 2 ,r',e). Be- 
sides, we keep only the contribution to the self-energy 
with one propagator of the superconducting fluctuations 
as illustrated in Fig. (3) 

£<->(r,r', £ ) = - 4 J gG>'<(r',r, W - e )^>(r,r', W ); 

(39) 

E fl ^(r,r',e) = -i J ^ [G< (r', r, u - e)L R > A (r, r', uj) 
+G A > R (r',r,u J ~e)L<(r 1 r' 1 u)]. 



The propagator of the superconducting fluctuations (see 
the end of Scc.[TTJ is determined by the standard geomet- 
rical series L(ui) = [— A -1 + IL(ui)]~ 1 , where II is approx- 
imated by the particle-particle polarization operator as 



CO 




co-e 



shown in Fig. [4] 

n< > (r,r',e) = - l J ^G<->(r, r', w - e)G<>> (r, r', e); 

(40) 

n^(r,r',e) = -i J [G<(r, r', u - e)G R > A (r, r', e) 

+G R ' A (r,r',w-e)G > (r,r',e)] . 

One may check that at equilibrium H 1 ^ = (1 + 
2np(uj))(U R q — U A q ), where np(uj) is the Bose distribu- 
tion function. After averaging over the disorder, H R and 
H A g are given by the standard expressions. 

We may now obtain the leading order corrections in 
the interaction to the electric current as a response to 
a temperature gradient in the linear regime. We should 
consider all possibilities to linearize the expressions for 
S V t and n V T with respect to VT in Eq. The 
diagrammatic interpretation for the different contribu- 
tions to the transverse electric current obtained in the 
quantum kinetic equation technique corresponds to the 
three diagrams shown in Fig. [5] After averaging over the 
disorder the leading contributions to the Nernst signal 
in the diffusive regime are obtained from the diagrams 
with three Cooperons^ presented in Figs. [Ha) and[HKb) 
and the Aslamazov-Larkin diagram^ shown in Fig. EUc). 
[The Cooperon is a singular diffusion propagator which 
describes the rescattering on impurities in the particle- 
particle channel.] Since we generate these terms using 
the quantum kinetic equation, the analytic structure of 




FIG. 3: The self-energy in the first order with respect to the 

propagator of superconducting fluctuations before averaging FIG. 4: The geometrical series describing the fluctuations prop- 
over the disorder. agator in the Cooper channel. 
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(a) 



(b) 



(cl 



FIG. 5: The diagrammatic contributions to the transverse com- 
ponent of je°" before averaging over the disorder. (The obvious 
counterpart diagram for (a) is not shown.) 




(VWWWWV. 



WlAAAOAAAT 



In the limit of small magnetic fields the Cooperons, 
like the quasi-particle Green functions, can be sepa- 
rated into the phase exp{2ie f*, A(r 1 )dr 1 /c} and the 
gauge invariant part at H = 0, C R ' A (p, e, u> — e) = 

[=Fi(2e - uj)t - DV 2 p ] ~\ see Appendix C in Ref. [H At 
a finite magnetic field, one may express the gauge invari- 
ant part of the Cooperon propagator using the Landau 
level quantization: 



C% A (e, tu-e) = [=R(2e - lu)t + n c r(N + 1/2)}' 



(41) 



where N is the index of the Landau level. Similarly, the 
propagator of the superconducting fluctuations written 
in terms of the Landau levels becomes: 



(a) 



(b) 



Cc) 



FIG. 6: The diagrammatic contributions to the transverse com- 
ponent of the Diagrams (a) and (b) describe the fluctua- 
tion of the superconducting order parameter decorated by three 
Cooperons and (c) is the Aslamazov-Larkin diagram. (The ob- 
vious counterpart diagrams for (a) and (b) are not shown.) 



n M = — 



In — +ip R<A (w,N) 



i) R:A (w, N) = if> 



•5 



i 

2 T 



ioj n c (N+i/2y 



AttT 



4ttT 



<rw ; 

(42a) 
(42b) 



the diagrams is given by the equation. 

To get the explicit expression for the current we re- 
turn to the gauge invariant equilibrium Green functions 
g given in Eq. 1161 Since we restrict our calculation to the 
limit lo c t <C 1 (where lo c = eH/m*c is the cyclotron fre- 
quency of the quasi-particlcs) , we may neglect the depen- 
dence of g eq on the magnetic field entering through the 
Landau quantization of the quasi-particles states. There- 
fore, the entire dependence of the quasi-particle Green 
functions on the magnetic field is through the phase. 
Unlike the quasi-particlcs. the Landau quantization of 
the collective modes (both the Cooperons and the fluc- 
tuations of the superconducting order parameter) can- 
not be neglected because the quantization condition for 
these modes is £I c /Tq > 1, where in the diffusive regime 
S! c = 4eHD/c is the cyclotron frequency in the Cooper 
channel. Note that il c cx lu c (sft) 3> lu c because the 
product of the Fermi energy and the mean free time is 
assumed to be a large parameter. [In Q c the effective 
charge is equal to 2e and the diffusion coefficient D re- 
places l/2m because in the Cooperons and the fluctua- 
tions propagators the term Dq 2 substitutes the kinetic 
energy p 2 /2m.] 



Here, i/j(x) is the digamma function. The primary goal 
of this calculation is to analyze the measurements of the 
Nernst effect in superconducting films££ In such films 
the electron states are not quantized and therefore v is 
the density of states of three dimensional electrons (as 
well as D). The parameter <; oc 1/(i/\sf) is important 
for understanding the difference in magnitude between 
the longitudinal and transverse Peltier coefficients. The 
longitudinal Peltier coefficient, a xx , contains an integral 
over the frequency that vanishes when c; = while the 
integrand determining a xy remains finite even in the ab- 
sence of As a result, in the expression for the Nernst 
coefficient given in Eq. [2] the second term in the numer- 
ator is smaller than the first one by a factor of the order 
T/(p\e F )M 

Using the expressions for the quasi-particle Green func- 
tions, the Cooperons and the propagators of the super- 
conducting fluctuations in the equilibrium state we may 
investigate the contributions of Gvt and Lvt to the 
current. Recall that we are interested in the transverse 
current. For illustration, let us show how to find the 
transverse current for one representative term out of the 
few contributions to the Aslamazov-Larkin diagram: 



ifS(rO = lim + ^ - + lim (Yl _ _ Yl _ 

2Tq J (2ir) 6 J ri2^ri\^2m 4mc 2m Amc J r 6 ^r 7 y 2m 4mc 2m Amc J 

go(*u r 2> ^SoVu, r 2 , w - e)g R (r lu r 12 , e)C R (r 2 , r 3 , e, u) - e)C R (r 10 , r U) e, u) - e)L R (r 3 , r 4 , w)Lf 9 (r 9 , r 10 , u) 

g R (r 5 ,r 6 ,e')g^(r 5 ,r s , e')g R {r 7 , r 8 ,e')C R (r 4 , r„ e', w - e')C R {r s , p„, e', uj - e')F(e, e', u). (43) 
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In Fig. [7] we indicate the spatial coordinates corresponding to the expression given above. Since in this part of the 
calculation we concentrate on the integration over the spatial coordinates, we collect all the frequency dependent fac- 
tors into the function F(e, e', 10) = e [tanh (e/2T) — tanh ((e — u>)/2T)] tanh ((w — e')/2T) dnp{u>) / 'dcu and leave them 
aside for a while. 

Next, we rewrite the Cooperons and the propagators of the superconducting fluctuations using the basis of the 
Landau levels states, </3jv,n(r) = RN,n(r)e m ^ / \/2tt (where i?jv,n(^) are the generalized Laguerre polynomials). In 
addition, we separate the quasi-particles Green functions into the phases and the gauge invariant Green functions. 
Then, following the flux technique introduced in Ref . |2?1 we rearrange Eq. 1431 as: 

■con, \ _ eV jy T f dede'dlV ST^ f , , - i cH(r 11 -r 1 )x(r 1 ~r 2 )/2c -ieH(r 5 -r 6 ) X (r 6 -r s )/2c (aa\ 

3eAr,)- wl ] (27r)3 2^J dr 2 ...dr 12 e (44) 

/V* ieH{ yi - jfa) V x 12 ieH{ Vll - y 12 ) \ _ R A R 

I™ — + 9o (r x - r 2 ,e)g ( ril - r 2 , u) - e )g£ r u - r 12 ,e) 

ri2-»ri \2m 4mc 2m 4mc J 

a (5 - - S + i£2 fe £i ) - - - r -'» 

e -ieH(r 8 -r II )x(r II -r 2 )/c-ieH(r 2 -r 6 )x(r 5 -r 8 )/c^ jV!o(l . 2 _ rj^fa _ ril )C«(e, W - e)C&(e, W - e) 

L« (o;)L^(o;)C«(e', w - e')^(e', w - e')F(e, e', w ), 



where Ih = y c/2eH is the magnetic length in the 
Cooper channel. In the last step we used the orthogonal- 
ity of the generalized Laguerre polynomials (an example 
for the treatment of the propagators in this basis can 
be found in Ref. l29l ). The first two exponents in Eq. [44] 
contain the magnetic fluxes accumulated in the triangles 
(r 1 ,r 2! i , i 1 ) and (r 5 ,r 6 ,r 8 ), respectively. One way to get 
the transverse current, is to extract the magnetic field 
from these two fluxes or from the diamagnetic terms. As 
a result the transverse current appears with the coeffi- 
cient uj c t. We neglect these terms; we will see that when 
the magnetic field responsible for turning the current to 
the transverse direction is extracted from the Cooperons 
or the propagators of the superconducting fluctuations 
one gets a much larger factor of the order Clc/T. There- 
fore, the integration over the coordinates of the two trian- 
gles can be done with the quasi-particle Green functions 
taken at H = 0: 



contribution to the current becomes: 



■con 
Je x 



eV y T 



8tt 2 T, 



'I" 

ox 



S-fV f dede'duo [ dr ^ (45) 
d ieHx\ , . 

— J <PM,0[ T ) 

A 1, ,\riR, 



ieHy 



<PN,o( T ) 



2D [ — 

ay 



C*(e',tJ-e')F(e,e',u;). 



The integral over the coordinate corresponds to the ma- 
trix element of the velocity operators (N, 0\V x V y \M, 0), 
where \M, 0) = ipM,o is the quantum state of a parti- 
cle with a mass equal to 1/2D in the M Landau level 
and zero angular momentum in the z-direction. Using 
the known properties of the Laguerre polynomials, the 
matrix element can be written as (N, 0\V x V y \M, 0) = 
2iD 2 [{N + l)5 NM -i - (M + 1)6 m ,n-i}/£ 2 h . Finally, the 



j e -2(r0 = -i eV « T 



^L-i^DV / dede'dw J2 ( N + !) 



C§(e, u> - e)C§ +1 (e, w - e)C§ (e> - e')C* +1 (e' ; u - e') 
[L%(u>)L% +1 {u) - L% +1 (u>)Lfr( U )] F(e, e',u). (46) 

In the limit H — > when the quantization of the collec- 
tive modes can be neglected, one may replace the Cooper- 
ons and the propagators of the superconducting fluctua- 
tions in Eq. 23] by the product of the phase terms (with 
charge 2e) and the corresponding propagators in the ab- 
sence of a magnetic field. Then, the contribution to the 
current at vanishingly small magnetic field can be found 
by employing the flux technique of Ref. [27l . One may 
check that the same result is obtained when the trans- 
formation from the discrete sum into an integral over a 
continuous variable is performed in Eq. |4"1)1 

Let us conclude with a remark regarding the diagram- 
matic interpretation of the different contributions to 
As already mentioned, the analytical structure and the 
expressions for the vertices of these diagrams were found 
from the quantum kinetic equation. In principle, the 
same diagrams can be calculated using the Kubo for- 
mula. However, if for simplicity one uses in the Kubo 





u ^\aaaaaaaaaJ 

r io 



FIG. 7: The Aslamazov-Larkin diagram. 
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formula the heat current operator of non-interacting elec- 
trons described in Eq. 0J the resulting expressions for 
these diagrams differ from those obtained in the quan- 
tum kinetic approach. Most important, as one can see 
from Eq. 1381 in the quantum kinetic approach the fre- 
quency accompanies the renormalized velocity, so that 
the expression for the electric current is generally of the 
form eg(e)vi(e)g(e)eVj(e)'VjT/T. In other words, the 
frequency appears together with the velocity that was 
already renormalized by the interaction. On the other 
hand, owing to the fact that the frequency in the sim- 
plified version of the Kubo formula is attached to the 
external vertex before the renormalization of the veloc- 
ity, the expression for the current has a totaly different 
structure. 

This is also the proper place to explain what is so 
unique in the superconducting fluctuations in the diffu- 
sive limit that leads to the giant Nernst effect. As any 
transverse current coefficient, a xy contains a difference 
of two almost equal terms. In addition, like all thermo- 
electric coefficients the integral over the frequency in a 
contains a factor of the quasi-particle frequency. Con- 
sequently, as discussed in Appendix [C] the contribution 
of the quasi-particles to the transverse Peltier coefficient 
includes two small parameters. The first, is the usual 
lu c t that appears in all transverse currents. The second 
is a reminiscence of the fact that the frequency factor 
(that in the Boltzmann equation is converted into the 
energy) is responsible for the vanishing of the Peltier co- 
efficient under the approximation of a constant density 
of states. When a non-constant density of states is con- 
sidered, the integration over the energy yields another 
small parameter proportional to T/ep. Now we turn to 
the contribution of the superconducting fluctuations to 
the transverse component of jg°™, and consider Eq . |4"B1 as 
a representative example. For the moment we ignore the 
factor e in F(e,e',u>) associated with the thermoelectric 
current. Then, the difference between the two almost 
identical terms results in an odd integrand with respect 
to the frequency of the superconducting fluctuations, w, 
which potentially may lead to the vanishing of a xy . So, 
how can the superconducting fluctuations induce a strong 
Nernst signal? The explanation lies in the fact that the 
Coopcrons accompanying the superconducting fluctua- 



tion depend on the frequency of the incoming/outgoig 
quasi-particles and not only on the frequency to carried 
by the fluctuations (Eq. HTj) . The dependence of the 
Cooperons on e combined with the frequency factor e 
in F(e,e',u)) save the situation. This is because the in- 
tegration over e results in an integrand that is an even 
function of to and, hence, there is no longer danger that 
the transverse Peltier coefficient vanishes. We shall see 
that instead of the two small parameters obtained for 
the quasi-particles, the contribution of the superconduct- 
ing fluctuations includes only one. Because of the extra- 
sensitivity of these fluctuations to the magnetic field this 
parameter is VL C /T$. 

V. FINAL EXPRESSIONS FOR THE 
TRANSVERSE COMPONENT OF \ c ° n 

If one examines Eq. [35] which presents the general ex- 
pression for the contributions to the electric current from 
Gvr and Lvt, one may notice that not all the terms 
contain the derivative of a Fermi distribution function. 
As one may expect, the terms in which the Fermi dis- 
tribution function is not differentiated contribute only to 
the transverse component of ] c e on and not to the longi- 
tudinal one. After integration over the Fermion degrees 
of freedom (the frequency e and the coordinates of the 
quasi-particles Green functions), the terms proportional 
to dnp(e)/de give two non- vanishing contributions. The 
first one corresponds to the Aslamazov-Larkin diagram 
presented in Fig.[SJc): 

[4( w )4 +1 ( w )-i« +1 ( w )4(^)] 

[V* (w, N) - i> R (to, N + I) + ip A (to, N) - (w, N + 1)] 

[Q C (N + 1/2) (V>« (w, N) - i) A (to, AO) 

-Q C (N + 3/2) & R (w, N + 1) - if> A (u, N + 1))] , (47) 

where the upper index in j conl enumerates the contribu- 
tion to the current, and Sij is the anti-symmetric tensor. 
The second contribution generated by terms with the 
derivative dnp(e)/de corresponds to the diagram with 
three Cooperons shown in Fig. [BJa): 



■con2 

Je i c 



eV,T 



lJ 47T 2 T 

iuj + Q c (N + 3/2) 



dw^n c {N + l)l-L*(u;) 



N=0 



1 R . dn P (u>) 



iu + fl c (N + 1/2) 
4^ 



(48) 



--L N n P (uj) 



ip A (w, N) - ip R (w, N) - ip A (w, N + 1) + V* (w, N + 1) 



icu + n c (N + 1/2) 
(4ttTo) 2 



ioj + n c (N + 3/2) 



47rT f! c 



+ N <-> N + 1 } +c.c 



r 



Here ip' R A an d V'h a correspond to the first and second The notation N <-> N + 1 means that N is replaced by 
derivatives of the digamma function defined in Eq. I42bl 
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N + 1 and the other way around in all the terms inside 
the curly brackets. Notice that there are no contribu- 
tions proportional to the derivative of the distribution 
function which can be attributed to the diagram shown 
in Fig Eft). 

Next, we discuss the group of terms that are pro- 
portional to np(e). The diagrammatic interpretation of 
these terms, which are generated by Eq. 1381 includes all 
three diagrams presented in Fig. El However, one may 
check that the contributions from the diagrams shown in 
Fig. Eft) andEft) are cancelled by a part of the contribu- 
tion from the diagram given in Fig. EI a ) ■ The remaining 
contribution is: 



7 -con3__ fe ..eV J -T 
Je i — tc 



-v / did (N + l)up(w) 

J N=0 



(49) 



v 4^ 2 T C 

{ [L*(uj) + L§ +1 (uj)} [i> R (to, N) - r (w, N + 1)] 
+^i*(o,)^ ( W> N) + ( W> JV + 1) J 

In the derivation of the different contributions to j™" we 
used the following identities for products of the distribu- 
tion functions: 



n F {e)nF{ui — e) = np(u>)[n F (e — u>) — nj?(e)]; (50) 

w) - n_p(e)] 

-np(u). 



dn F (uj — e) , . dn P (uj) 

-n F (e) = — = [n F (e 



<9w 



9lj 

dn F (oj — e) 



div 



Further analysis of j^ on at arbitrary temperatures and 
magnetic fields shows that in j^"? 1 and j™" 2 the integra- 
tion over the frequency accumulates at u ~ T <C 1/r. 
As a consequence of the narrow range of the integration, 
the final expressions for these two contributions vanish 
in the limit T — > 0. In contrast, in j™™ 3 the integration 
over the frequency is not limited to small frequencies and, 
hence, the outcome of the integration depends logarith- 
mically on the scattering rate 1/r which acts as an ul- 
traviolet cutoff. In addition, as the temperature goes to 
zero there is even a more serious problem with this term, 
because its pre- factor is proportional to Q c /Tq. Such a 
dependence on the temperature violates the third law of 
thermodynamics. [The connection between the third law 
of thermodynamics and the Nernst effect was discussed in 
the Introduction.] We shall see that the dangerous parts 
in jg°f 3 are cancelled out by the magnetization current 
that up to now we have not yet considered. 



VI. THE MAGNETIZATION CURRENT AND 
THE THIRD LAW OF THERMODYNAMICS 



that this part of the Green function does not contribute 
to the magnetization current. On the other hand, the 
explicit dependence of the local equilibrium Green func- 
tion on the center of mass coordinate leads to a non-zero 
contribution to the magnetization current: 



■mag 
Je 



RVT 

Tn 



(51) 



dc 



Thus, Gvt and Gi oc ~ eq are complementary to each 
other; while the first contributes only to jg ™, the other 
one fully determines j" lag . One should recall that we are 
looking for a current that does not vanish after spatial av- 
eraging, i.e., after integration with respect to the center 
of mass coordinate R. Since in the process of averag- 
ing over R we may integrate by parts, the magnetization 
current can be written as 



■mag = 2 i £ljC M z lim 



2tt T 



9fq(P,£\ A, imp). (52) 



Here we integrated by parts over the frequency as well. 
One may recognize that j™ a f is directly related to the 
magnetization density at equilibrium: 



To ' 



(53) 



The result demonstrates the strength of the quantum ki- 
netic approach. This method provides a way to derive 
both components of the current without engaging any 
thermodynamical arguments. 

Actually, at this point one may employ in Eq. [53] the 
known expression for the magnetization in the presence 
of superconducting fluctuation. Still, since we are inter- 
ested in the interplay between the quasi-particle excita- 
tions and the fluctuations of the superconducting order 
parameter, let us derive the expression for the first order 
correction to the magnetization induced by the fluctua- 
tions starting with Eq. [52j Using the standard identities 
for the Keldysh Green function at equilibrium, one gets: 



- 7 ^ lim 

T r'^r 



de 
2^ 



dr t dv 2 [M z (v) + M z (r')] 



(54) 



x n F (e) [go (r - r 1; e; A, imp)cr A (r 1 - r 2 ,e; A, imp) 
xgoVa - r', e; A, imp) - gg(r - r lt e; A, imp) 
xa R (r 1 - r 2 ,e; A, imp)g^(v 2 - r' , e; A, imp)] . 



In this section we examine the magnetization current 
given in Eq. [37l In general, we need to insert the VT- 
dependent part of the Green function, Gi oc ^ eq + Gvt 7 
into Eq. 1371 Since after the averaging over the disor- 
der Gvr(r — > r', e) is translationally invariant, it is clear 



Here, for convenience, we returned to the initial coordi- 
nates. Next, we use the fact that the equilibrium Green 
function in the absence of fluctuations satisfies the follow- 
ing identity — dg R A /dH = g^" A 'M.g R ' A . Therefore, the 
expression for the magnetization current can be rewritten 
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FIG. 8: An illustration of the relation between the magnetiza- 
tion current term that is obtained from the local equilibrium 
Green function and the thermodynamic diagram for the mag- 
netization. 




as 



■ mag 

3 e i 



-lice 



: To 



lim 



2^ 



dr 1 np(e) 



(55) 



d9o(r- r 1; e; A, imp) A 



dH z 



a (ri — r', e; A, imp) 



dg^{r-r u e;A,imp) R 



dH z 



a (r x - r 2 , e; A, imp) 



Finally, using the explicit expression for the self-energy 
and rearranging all the terms we reformulate the expres- 
sion for the magnetization current in terms of the prop- 
agator of the superconducting fluctuations: 



■mag 

3e i 



-ice in — ; — lim 

1 r'^r 



2^ 



dr 1 np(w) 



(56) 



dH R (r — r l5 ui; A, imp) 
dH z 

dH A (r — t-i,u)\ A, imp) 



-ice. 



dH z 
V T d 



lim 



L' R (r 1 — r', w; A, imp) 
L A (r x - r 2 ,w; A, imp) 
-np(w) 



[lnL^ 1 (p, u>; A, imp) — lnZ/^ 1 (p, w; A, imp)] . 

The transition from Eq. to the last line in Eq. [5(5] is 
illustrated in Fig. [H Averaging over the disorder and 
transforming from the expression for the propagator as a 
function of the coordinates to the basis of Landau levels, 
one obtains the known expression for the correction to 
the magnetization in the lowest order with respect to the 
fluctuations^^ 



■maq 
Je i = E 



~V J 2^i ^ Up ^ 

J N=0 



1 T dH 



(57) 



The discussion of higher order corrections to the magne- 
tization is given in Ref. l30l . 

Similar to j™" 3 in Eq. 221 the integration over the fre- 
quency in the magnetization current is restricted by the 
scattering rate, and at low temperature j™ ag also diverges 
as fl c /T. The opposite sign of the magnetization current 
relative to j™' 13 suggests that these dangerous parts may 
cancel each other making the Nernst signal compatible 



FIG. 9: The phase diagram for the Peltier coefficient a xy . We 
indicate the equations in the text which give the corresponding 
expressions for a xy in the different limits. f2 c = 4eHD/c is the 
cyclotron frequency for the fluctuations of the superconducting 
order parameter in the diffusive regime. 



with the third law of thermodynamics. Another hint for 
this cancellation is the similar analytical structure of j™" 9 
and j™" 3 . [All the terms in Eq. g51 and in Eq. \E7\ are a 
product of either retarded or advanced functions only] 
In Appendix [B] we show that the diverging parts of the 
magnetization current indeed identically cancel out the 
diverging parts of j™" 3 . We demonstrate that the to- 
tal current is independent of r in the whole temperature 
range T <C 1/r. As a result, the Nernst signal is regular 
at T — > 0. Moreover, the contributions which are con- 
stant with respect to the temperature also vanish and the 
remaining terms are linear in T. 



VII. THE PHASE DIAGRAM FOR THE 
NERNST EFFECT - COMPARISON BETWEEN 
THE THEORETICAL RESULTS AND THE 
EXPERIMENT 

In the following part we present the theoretical expres- 
sions for the transverse Peltier coefficient for a supercon- 
ducting film in the normal state for various regions of 
the temperature and the magnetic field. The phase dia- 
gram for the Peltier coefficient is plotted in Fig.[5J In the 
area below the line ln(T/T c (ff)) = Sl c /4.itT the Landau 
level quantization of the superconducting fluctuations be- 
comes essential. The line \n(H/H C2 (T)) = 4irT/fl c sep- 
arates the regions of classical and quantum fluctuations. 
From now on Tq is replaced by T which represents the 
spatially-averaged temperature. 

For a small magnetic field, il c <C T, close to the tran- 
sition temperature (T ss T c ) the leading contribution to 
a xy is given by the Aslamazov-Larkin term (see Fig. [6] 
(c)) and the magnetization current: 



*xy 



192T In T/T C {H) 



(58) 



In the previous section we discussed in details the impor- 
tance of the magnetization current in cancelling the quan- 
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turn contributions to the Nernst signal. In the vicinity 
of T c one can interpret the expression in Eq. [55] in terms 
of the classical picture in which the Cooper pairs with a 
finite lifetime are responsible for the thermoelectric cur- 
rent. The magnetization current is just equal to —2/3 
of the leading order contribution from the Aslamazov- 
Larkin term. Note that unlike the electric conductivity, 
a xx , for which the anomalous Maki-ThompsonA£ and the 
Aslamazov-Larkin terms yield comparable corrections, 
the contribution from the anomalous Maki-Thompson 
term to the Nernst signal is ~ (T/ep) 2 -C 1 smaller 
than the one given by Eq. 1581 Therefore, it is natural 
that in the vicinity of T c our result coincides with the ex- 
pressio n 13 ' 31 obtained phenomenologically from the time 
dependent Ginzburg-Landau equation (TDGL). 

When temperature is increased further away from the 
critical temperature, the sum of the contributions to the 
transverse Peltier coefficient from all the diagrams and 
the magnetization current yields: 



ixy 



247r 2 TlnT/T c ' 



(59) 



A comparison between the transverse Peltier coefficient 
in the vicinity of T c (Eq. I58[) and far from the transition 
(Eq. IB"9"|) reveals that the two expressions differ only by a 
numerical coefficient. The similarity between the expres- 
sions for a xy in the two different limits is not seen in para- 
conductivity. This is a consequence of the cancelation of 
the quantum contributions to by the magnetization 
current, which is specific for the Nernst effect. Away 
from the critical region r>T c , the quantum nature of 
the fluctuations reveals itself in contributions to j™" 3 and 
ynag fag^ CO ntain an integration over a wide interval of 
frequencies between T and 1/r. As a result, these terms 
become of the order In (In 1/TY) — In (In T/Tc). However, 
as we show in Appendix [B] these r-dependent terms in 
jcon3 anc j ynag cance i eacn other <22 The Peltier coefficient 
far from T c demonstrates how the third law of thermody- 
namics constrains the magnitude of the Nernst signal not 
only at T — > but also at high temperatures, T ^> T c . 

The comparison of our result with the experimen- 
tal observation of Ref. [H for two Nbo.15Sio.85 films of 
thicknesses 35nm and 12.5nm (with critical tempera- 
tures T c = 380mK and T c = 165mK, correspondingly) 
is given in Fig. [TU1 The Peltier coefficient depends on 
the mean field temperature of the superconducting tran- 
sition, T c , and on the diffusion coefficient through O c . 
Throughout the paper we fit the data using the same 
diffusion coefficient D = 0.187cm 2 /sec which is within 
the measurement accuracy of the value that is extracted 
from the experiment in Ref. 0. 

The cancelation of the terms proportional to O c /T in 
the limit T <C fi c was described in the previous section. 
(Without this cancellation we would get a finite Nernst 
effect in the limit T — > and the third law of thermo- 
dynamics would be violated.) After the cancellation, the 
remaining contributions to a xy in the limit T — > are 
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FIG. 10: The transverse Peltier coefficient a xy divided by the 
magnetic field If as a function of In T/T c in the limit H — ► for 
films of thicknesses (a) 35nm and (b) 12.5nm. The experimental 
data of Ref. is presented by the black squares and the solid 
line corresponds to the theoretical curve given by Eq. 1591 The 
insert presents the fitting of the data in the vicinity of T c with 
Eq.EH 



linear in the temperature: 
eTln3 



~3n c (lnH/H C2 {T))< 



for H w H, 



c 2 • 



and 



2eT 



a. 



3fl c In H/H c . 



for if > H c 



(60) 



(61) 



Similar to the limit Q c < T the integrals determining 
the final expression for a accumulate at low frequency. 
This situation is not typical for fluctuations induced by 
a quantum phase transition. Notice that a xy changes 
its sign in this region. Since the transverse signal is 
non-dissipative the sign of the effect is not fixed. As 
mentioned before, in the vicinity of T c for Q c <C T, 
the main contribution to the Peltier coefficient is from 
the Aslamazov-Larkin term and the magnetization cur- 
rent. The magnetization current is opposite in sign to 
the Aslamazov-Larkin terms and equals 2/3 of it. When 
crossing to the region \a{T /T C {H)) < Sl c /T (see the 
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phase diagram in Fig. [9]) the contribution from the mag- 
netization current grows. To the first order in £l c /T 
the magnetization current cancels the Aslamazov-Larkin 
term and the Peltier coefficient turns out to be propor- 
tional to 0[(il c /T) 2 ]. Lowering further the temperature 
and increasing the magnetic field one reaches the re- 
gion Q c > T and ln(H/H C2 (T)) < T/tt c . In this region 
the magnetization current becomes dominant. Since the 
magnetization current gives a contribution that is oppo- 
site in sign to the Aslamazov-Larkin term, we obtain that 
the Peltier coefficient is negative. 

In Fig. QT]we plot the Peltier coefficient for the 35nm 
film as a function of the magnetic field at a temperature 
higher than T c . We take T C MF = 385mK to be slightly 
higher than the measured critical temperature anticipat- 
ing a small suppression of the temperature of the transi- 
tion by fluctuations. [The data in Fig.QT] unlike the data 
in Fig. 10(a) is presented in linear rather than a logarith- 
mic scale. Therefore, this fit is much more sensitive to 



x 10" 



the input parameters compared to the one in Fig. 10(a) 
As a result, the small deviation of T^ F from the mea 
sured T c can be noticed. For consistency, we use the same 
value of T^ F also in Fig. 10(a) ] Fig. [TT1 demonstrates 



the agreement between the theoretical expressions and 
the experimental observation for a broad range of mag- 
netic fields. In addition, we show that the experimental 
data is well described by Eq. [55J in the limit of vanish- 
ing magnetic field; see inset of Fig. [TTJ Since Eq. [55J is 
valid in the limit il c <C T, it can describe only the first 
few point in the measurement. In order to fit the en- 
tire range of the magnetic field we had to include higher 
order terms in Cl c /T. For that we needed to sum the 
contributions from all diagrams and the magnetization 
current. We performed the calculation assuming that 
ln(T/T c (H)) <C 1; therefore the theoretical curve starts 
to deviate from the measured data when ln(T /T C (H)) is 
no longer small (H f=a lTesla). 



VIII. SUMMARY 

We demonstrated that the contribution from the fluc- 
tuations of the superconducting order parameter to the 
Nernst effect in disordered films is dominant and can 
be observed far away from the transition. We showed 
that the important role of the magnetization current is 
in cancelling the quantum contributions, thus making the 
Nernst signal compatible with the third law of thermo- 
dynamics. The third law of thermodynamics constrains 
the magnitude of the Nernst signal not only at low tem- 
peratures, but also far from T c . As a consequence of this 
constraint the phase diagram is less rich and diverse than 
one expects in the vicinity of a quantum phase transition. 

The Nernst effect provides an excellent opportunity 
to test the use of the quantum kinetic approach in the 
description of thermoelectric transport phenomena. We 
showed that in this scheme we get automatically all con- 
tributions to the Nernst coefficient as response to the 




H (Tesla) 



FIG. 11: The transverse Peltier coefficient a xy as a function of 
the magnetic field measured at T = 410mK. The black squares 
correspond to the experimental data of Ref. 8 while the solid 
line describes the theoretical result. The arrow on the phase 
diagram illustrates the direction of the measurement. In the 
insert the low magnetic field data is fitted with the theoretical 
curve given by Eq. [58] 



temperature gradient, in particular the one from the 
magnetization current. This is an advantage of the quan- 
tum kinetic approach but it is not the only one. This 
method also allows us to verify the Onsager relations; 
a comprehensive discussion of this issue is presented in 
Appendix [X] The fact that the we were able to find in- 
dependently the two off-diagonal components of the con- 
ductivity tensor, <xy and a^j, and verify that they are 
connected through the Onsager relation assures that the 
quantum kinetic approach developed in this paper and in 
Ref. [l3 gives the correct expressions for the electric and 
heat currents. 

Finally, we should remark that our results for the 
Peltier coefficient differ in few aspects from those ob- 
tained recently in Ref. [33| using the Kubo formula. As 
we already discussed in the end of Sec. IIV1 the simplified 
Kubo formula cannot give the correct electric current as a 
response to a temperature gradient. Therefore, the claim 
of the authors of Ref. |33| that the difference between the 
Nernst signal calculated using the simplified Kubo for- 
mula and the quantum kinetic approach is only in the 
numerical coefficients is unacceptable. The expression 
given in Ref. [33] for a xy in the vicinity of T c cannot fit 
the experimental data, and it also contradicts the phe- 
nomenological result of the TDGLi^i The only fit of the 
experimental data presented in Ref. [HI is a logarithmic 
plot of the Nernst signal as a function of temperature 
using the formula for temperatures not too close to T c . 
Such a logarithmic plot is not very sensitive to the nu- 
merical coefficients. The striking agreement between our 
results and the experimental data, in particular our abil- 
ity to obtain the non-trivial dependence of the Nernst 
signal on the magnetic field and the fact that we repro- 
duced the phcnomcnological resuh^I reinforces us in the 
correctness of our method. 
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APPENDIX A: ONSAGER'S RELATIONS 

In this Appendix we compare the electric current aris- 
ing as a response to a temperature gradient and the heat 
current generated by an electric field. We verify that for 
the Gaussian fluctuations of the superconducting order 
parameter, the two expression are connected through the 
Onsager relations^ &ij(B) = Xb<Xji(— B). [In a similar 
way, in Ref. [13 we demonstrated the Onsager relations for 
the longitudinal current in the presence of the Coulomb 
interaction.] 

The derivation of the electric current induced by the 
temperature gradient was presented in Sec. [ill where we 
discussed the two contributions to the current. The first 
one, jg°™, was found using the continuity equation and its 
expression before the expansion in the superconducting 
fluctuations is given in Eq. 1381 The second contribution, 
analyzed in Sec. IVI1 is from the magnetization current. 
This term contributes only to the transverse current, and 
it can be written as j™" 9 = ce itj M z (-V 'jT '/T ). 

Next, we sketch the derivation of the heat current as a 
response to a uniform electric field. For this purpose we 
use the heat continuity equation, Q(r,t) + Vj£ on (r, t) = 
jg°"E. The product jg OTl E describes the work performed 
by the electric field on the current. Since the electric field 
cannot perform any work on the magnetization current, 
only enters to the RHS of the continuity equation^ 
The heat density in the presence of an electro-magnetic 
field is a function of the magnetization and the electro- 
chemical potential: 

dQ(r, t) = dh(r, t) - (fx- etp(r))dn(r, t) + MdH(r, t), 

(Al) 

where h(r, t) is the Hamiltonian density. To find the time 
derivative of the magnetic field, we turn to the Maxwell 
equation H = — cV x E. Thus, the heat current is de- 
scribed by the following equation: 

Vj ft = -dh(r, t) + - e(p)dn(r, t) - 'f e on Vtp + cMV x E 

(A2) 



In Ref. [I?] we showed that for H = the expression for 
the heat current found from the continuity equation is 



•con 
Jh 



(H = 0) = lim {d t + ieip{r) - d t > + ie<p(r')) 



r — >r 

t'-*t 



This result was obtained for the Coulomb interaction 
but it also holds for the superconducting fluctuations 
because, unlike the charge, there is no principle differ- 
ence between the way the fluctuations in the density and 
Cooper channels carry heat. Here, we use the fact that 
in the presence of an interaction field, such as A, which 
does not have its own dynamics, the heat current can 
be formulated in terms of the quasi-particle Green func- 
tion alone. This is compatible with the observation that 
according to the kinetic equation given in Eq. 1231 the 
temperature gradient is coupled to L only through the 

quasi-particle Green functions inside II. [When interac- 
tions with dynamic fields like phonons are studied the 
heat current acquires additional terms.] 

Although we restrict our derivation to the regime of lin- 
ear response with respect to the electric field, the source 
term jJ; on E is still important. This is because the source 
term makes the heat current to be gauge invariant as 
Eq. I A3I reveals. In principle, there may be an additional 
contribution to the heat current from the charge current 
carried by the superconducting fluctuations (correspond- 
ing to the RHS of Eq. [29]) . This issue is not addressed 
here because such contributions are beyond the linear 
response. 

When we consider the effect of applying a magnetic 
field the expression for the heat current given in Eq. IA3I 
has to be modified. The first change is simply to in- 
clude the vector potential in the velocity as it is shown 
in Eq. [22j We denote the contribution from the heat cur- 
rent given in Eq. IA3I with the modified velocity as j™" 1 - 
Besides, there is an additional contribution to the heat 
current, j£ on2 , from the last term in Eq. lAll that contains 
the magnetization: 



cM(V x E) = cV(E x M). 



(A4) 



Here we used the fact that under the condition of a con- 
stant magnetic field V x M = 0. For the setup of the 
Nernst measurement (see Fig. [T]), in which the magnetic 
field is aligned along the z directions, the contribution of 
the magnetization to the heat current is 



(A5) 



dv 1 dt 1 Yiv.t^v^t^Giv^t^v' ,t') + h.c. (A3) 



At this stage, one may wonder whether there is a con- 
tribution to the transverse current that cannot be found 
from the continuity equation, i.e., a term of the form 
V x W. In the case of the electric current generated 
by the temperature gradient, we saw the term of this 
kind is the magnetization current. This term does not 
vanish and contributes to the transport electric current 
because the non-uniform temperature induces a coordi- 
nate dependent magnetization (see Sec. IVI[) . However, 
in the presence of a constant electric field, the system re- 
mains uniform. Therefore, the quantity that we denoted 
by W should be independent of the spatial coordinate 
and, hence, V x W = 0. 

Let us compare between the electric current as a re- 
sponse to a temperature gradient and the heat current 
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generated by a electric field when a magnetic field is ap- 
plied. One may immediately notice that j£°™ 2 and j™ a f 
(given in Eqs. IA5I and [53] respectively) satisfy the On- 
sager relations: 



■ mag 

le j 



(-H) 



-V,T 



(A6) 



It is interesting that these two terms coincide although 
they seem to have a different origin. The contribution of 
the magnetization, V x M, to the electric current cannot 
be found from the continuity equation, and it is non-zero 
due to the dependence of Gi oc -eq on the center of mass 
coordinate. When the response to a uniform electric field 
is considered, the Green functions arc independent of the 
center of mass coordinate. Still, there is an equivalent 
contribution to the heat current arising from the conti- 
nuity equation. 

Now, we show that j™" and j£ onl also satisfy the On- 
sager relations. In order to find j£ onl we need to know 
the expression for the electric field dependent propaga- 
tors. The electric field dependent Green function can be 



written in the following form>i£ 

G E (p, e; A, imp) = g eq (e) E E (e) g eq (e) 
ieE 



(A7) 



dgeg (e) - . ,^ , , dg eq (e) 

— Ve q {£)g eq (e) - geq (e) v eg (e)- 



dc 



ck 



The above equation is similar to Eq. [21] for Gvt ■ Since 
Eg contains also the electric field dependent propagator 
of the superconducting fluctuations, we have to find the 
equation for Le- Owing to the fact that the supercon- 
ducting fluctuations carry charge, their coupling to the 
electric field is more complicated than their dependence 
on the temperature gradient described in Eqs. 121)1 and |2"T1 



Le(p, A, imp) = —L eq (u) IIe (u>) L eq (u) (A8) 
+ zeE 



dL eq {u)~ f f - dL eq {uj) 

V eq {Uj)L eq (w) - L eq (u>) V e «j(w)- 



8lu 



du 



Inserting the expression for Ge into Eq. IA3I and extract- 
ing the lesser component, we get: 



Jut = ^ / ^^T 1 [ W f(e)4(e)<(e)4(e) + vf {e)g«{e)vf (e)4(e) - v*{e) g *{e)vf{e) g *{e) 



2tt de 



-gf q {e)vf{e)gf q {e)vt{e)]+eE j \ ^en F (e) 



de 



R/ \^5en( e ) Rf \ Rl \ Rr \ Rr \ Rl \ ^3en( e ) 

g £ <(e)<(e) - <(e)gf g (e)^(e) ^ 



* / ^^fW4( c ) [^(e)(l-n F (e)) + S>(e)n F (e)] (<£(e) - 4(e)) + cc. 



(A9) 



The fulfilment of Onsagcr relation demands micro- 
scopic reversibility, which in our case implies that 
G(r,r',e;H) = G(r',r,e;-H) and L(r,r',e;H) = 
L(r', r, e; — H). Since the currents contain a trace over 
the coordinates, it is obvious that the first two terms in 
Eq. IA9I and the first two terms in Eq. [35] arc connected 
through the Onsager relations. 

Let us examine the remaining terms in \ c ° n and j™" 1 - 
The electric current contains not only the contribution of 
Svt but also terms with LTvt, while the remaining part 
of the heat current contains only Se- Actually, when we 
treat £e we must consider the possibility that E enters 
also through L. Then, as we show later, since the equa- 
tion for Le includes more terms than the equation for 
-Lvt (sec Eqs. IA8I and [27]) . the corresponding contribu- 
tions to the electric and heat currents coincide. 

Let us start with the contributions to the currents in 
which £ depends on the electric filed/temperature gradi- 
ent through G rather than L. Since we are interested in 
the effect of Gaussian fluctuations, we can use the expres- 



sions for the Ge and Gvt in the absence of interactions: 



GUe) = Me)n F (e) 



g e v a9o \ e ) 9o V e ) v o g e 



+ \F 3 {e)^^g^eH [gfc) - g«(e)] - c.c, (A10) 



where Fj(e) is equal to eEj and eVjT/To, respectively. 
The only difference between the above equation for the 
lesser component of Gf and the one for the greater com- 
ponent is that in the latter the distribution function, 
rip(e) should be replaced by n F (e) — 1. (In Eq. lAlOl and 
below we start to place the spatial direction indices also 
as superscripts.) 

Using the identities given in Eq. I50[ one may notice 
that for the discussed contributions to the currents (de- 
noted as j e ,fc(Gp, H)) only the terms proportional to the 
derivative of the distribution function in Eq. IA10I do not 
vanish. As a result we get: 
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v,(r„, r,) 




v.fe.r,) 



FIG. 12: A diagrammatic representation of the contribution to 
the current written in Eq. IA11I For simplicity the scattering by 
impurities is not indicated. 



4(G Fl H)=-, 



deduj 



dr 2 ...dr 6 f e / h (e)Fj(u - e)^^,^) [g^r^r^e) - g A {r 1 ,r 2 , e)] 



(All) 



2 J (2tt) 2 

[L? q (r 2 ,r 5 ,u}) - Lf q (r 2 ,T 5 ,u)] ^(r 4 ,r 3 ) [ ff ^(r 5 ,r 4 ,cj - e) - g£ (r 5 , r 4 , u - e)] [g^(r 3 ,r 2 ,w - e) - 5( f (r 3 , r 2 , cj - e)] 
[.9?(i"6,r 1; e) -5^(r 6 ,r l5 e)] ^^[?i F (e - - n F (e)]. 



r 



Here, / e (e) = — e and /ft(e) = e. The diagram- <5(r — r'). Finally, we shall change the frequency as fol- 
matic representation of the above expression is presented lows e — > cj — e: 
in Fig. Q~2J Recall that the bare velocity Vg(r,r') oc 



3l/h(GF,n) = -\] j^ dr 2--dr 6 f e/h {uj - e)F J {e)v (r 6 ,r 1 ) [^(r^ r 2! a; - e) - ^(r 1; r 2 , w - e)] 
K(r 2 ,r 5 ,a;) - £^(r 2 , r 5 , w)] ug(r 4 ,r 3 ) [.g^(r 5 , r 4 , e) - 5 ^(r 5 , r 4 , e)] [g£(r 3 , r 2 , e) - .g^(r 3 , r 2 , e) 
[9o(r e ,T u u; - e) - gf(T e ,T u LJ - e)] [n F (e-w)- n F (e)] . 



(A12) 



Under the condition of a microscopic reversibility men- 
tioned previously, we see that the electric current created 
by a temperature gradient as given in Eq. IA11I (with 
f e (e) = -e and Fj(w - e) = (u - e)VjT/T ) and the 
heat current genereted by an electric field as described in 
Eq. IA12I (with fh(u) — e) = u — e and Fj(e) = Ej) satisfy 
the Onsagcr relations. Thus, the microscopic reversibility 
and the Onsager relations emerging from it correspond 
to reading the diagram in Fig. [T^] from right to left in- 
stead of left to right (i.e., reading it in Hebrew instead of 
English). 



Next, we shall examine the contribution to the cur- 
rents in which the propagator of the superconducting 
fluctuations (or the polarization operator) depends on 
the electric field/ temperature gradient. We start with 
the corresponding contribution to the heat current as a 
response to an electric field, j/ 1 (Le,H). Using the iden- 
tities for the distribution functions in Eq. 1501 this term 



can can be written as: 

f h (L E ,H) =j j^-^dr 2 ...dr 4 ev l (r 4 ,r 1 ) [n F (e) - n F (e - uj)} 

[deg Oi , r 2 , e) ~ 9e q ( r i > r 2 , <0] [s£ ( r 3, r 4) e) - g% (r 3 , r 4 , e)] 

[g£ q (r 3 ,r 2 ,uje) - g? q (r 3 ,r 2 ,u - e)] 

[(1 + n P (cj))i^(r 2 ,r 3 , uj) - u P (w)L^(r 2 ,r 3 ,w)] . 

(A13) 

Using the definition of II vt, one may notice that the 
expression for the heat current can be rewritten in the 
form: 



dw , / dnp{u) 
2^ \ dJ~ 



(A14) 



[ n v*T( r/ ! r > + np{u>)) ~ n v,r( r '' r ' w)n P (c 
(l+w P (w))ll(r,r» - n P {u)L% (r,r',w) 



Let us now turn to the equivalent contributions to 
the electric current created by a temperature gradient, 
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je(ivT,H). Performing manipulations similar to those in the expression for the heat current we get: 
I 

j£(L V T,H) = -ie J < ^-dv 2 ..dr l [(£*(r 3 ,r 4 ,u;) - L^(r 3 ,r 4 ,w))Vf (r i ,r 1 ,cj)L eq (r 1 ,r 2 ,u)) - L^(r 3 , r 4 , w)V/(r 4 , r l5 w) 

(A15) 

x(if g (ri,r 2 ,o;) - L^(r lf r 2 ,cj))] [(1 + n P (w))Ii< iT (r 2 , r 3 , cj) - n P (w)II> . T (r 2 , r 3 , w)] 

+ 6 / ^3^ dr 2" dr 4 u o( r *' r O [5^(i"i,r 2 ,e) - .g^(r 1; r 2 ,e)] [g^(r 3 ,r 3 , w - e) - g£(r a , r 2 , cj - e)] 

x [(1 + n P (u))L< , T (r 2 , r 3 , w) - n P {uj)L> , T (r 2 , r 3 , w)] [n P (e) - n P (e - w)] [sjf (r 3 , r 4 , e) - c^(r 3 , r 4 , e)] . 

I 

Keeping in mind the definition of the polarization op- ITe(w). Then we may collect the two contributions into 
erator, we may rewrite the second integral in terms of a more compact expression: 

I 



<9?ip(cj)\ 1 



dui 



1 f doj 

j^(i V T,H) = -— J — dr 2 ...dr i [(1 + n P (uj))U< iT (r 3 ,r 4 ,uj) - n P (uj)n> iT (r 3 , r 4 , u>)] 

- ieE i dnP d u [ L tg( r ^^)^( r ^2,co) (L*(r 2 ,r s ,u>) - L^(r a ,r s ,w)) - (if 9 (r 4 , r 1; w) - L^(r 4 , r 1; w)) 
Vf (r^r^Lf^r^r^o;)] + Lf q (r i ,r 1 ,u) (1 + n P (w))II|. (r^r,, u) - n P ((j)U>. (r x , r 2 , Xf 9 (r 2 , r 3 , w) I . 



(A16) 



The expression inside the curly brackets in the above 
equation can be written as (1 + np(uj))L^(uj) — 
np(w)Lg(w). To obtain this identity one should find 
the lesser and greater components of L from Eq. IA81 A 
simple calculation reveals that in this combination of L^. 
and only the terms proportional to the derivative of 
the Bose distribution function and those including hlg :> 
(which also contain dnp(ui)/duj) give a non-zero contri- 
bution. Once again, if we invert the direction of the 
propagation of all the ingredients in Eq. IA16I and also 
the direction of the magnetic field we get that the ex- 
pression for this last contribution to the electric current 
becomes: 



du) , ( dnp{ui) 



-4/1* 



V eu (A17) 

x [n^.. r (r', r, + n P {uj)) - n^ T (r', r, w)n P (w)] 
x (1 +n P (w))L^(r,r',u) - np(uj)L^. (r,r' ,uj) . 

Comparing j'^LviT, — B) given above and j\(Ej, B) pre- 
sented in Eq. IA141 one immediately sees that they are 
indeed connected by the Onsager relations. 

In conclusion, we demonstrated the Onsager relation 
for the Gaussian fluctuations of the superconducting or- 



der parameter (i.e., to the leading order in (e p t) _1 ). The 
structure of the expressions for the electric and heat cur- 
rents, Eqs. [33J and IA91 indicates that the same is true 
for any order. [An example for a general proof of the 
Onsager relations is given in Sec. 6 of Ref. hj.1 



APPENDIX B: CANCELATION OF THE 
THERMOELECTRIC CURRENT IN THE LIMIT 

T —> 

In this Appendix we demonstrate how the contribution 
to the transverse component of \ c ° n that does not vanish 
when T — > is cancelled by the magnetization current 
given in Eq. (53J From the three terms constituting jg 0n , 
which arc described in Sec. IIV1 the one that remains at 
low temperatures, j™" 3 , is presented in Eq.^SJ Let us re- 
store the general expression from which j™" 3 originates: 



■con3 _ eV J T f de 



To 



/da 
—dr 2 ...dr 6 n F (e)[v l (r 6 ,r 1 ) (Bl) 



x 9o(ri, r 2 , e)af q (r 2 ,r 3 , e)g^{r 3 , r 4 , e)^(r 4 , r 5 ) 3 ^(r 5 , r 6 , e) 
-«o(r 6) ri)go(T u r 2 , e)i^(r 2 , r 3 )g£(r 3 , r 4 , e)of g (r 4 , r 5 , e) 
Xfl^(r 5 ,r 6) e)] + c.c. 
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FIG. 13: (a) The contribution to j£ em (a) before and (b) after 
averaging over the disorder. 



Here we use the notation vo(r, r') for the bare veloc- 
ity (see Eq. 1221 with a eq taken to be zero). Following 
Refs. and l36l and recalling that we are interested in 
the current averaged over the space, we may replace 
the convolution g R (e)VQg R (e) in the above equation with 
— i(r — r')g R (r, r', e). Then, the contribution to the cur- 
rent from jg°™ 3 is: 



Obviously, when we add the magnetization current pre- 



sented in Eq.[52]to jg ™ 3 , only the second term in Eq. 
remains. We shall denote this remaining term as 'f e em . 

Next, we replace a eq with its explicit expression given 
in Eq. [33] 



■rem _ c f d ^ 

Je i ~ e 



<T f dedto , , 
— / , n v , dr 2 ...dr i np(Lu)E za p(r 2 + r 3 
o J (2tt)^ 



J 



con3 



. eVjT fde 

"^7— / —dr 2 ...dr 4 n F {e) 



(B2) 



w o( r 4> r i)( r i - r 2 )iS?(r 1; r 2 , e)crf g (r 2 , r 3 , e) ff ^(r 3 , r 4 , e) 

-t£(r 4! rOO*! - T 3 )jg^(T u T„ e)a R (r 2 ,r 3 , e)g^(r 3 , r 4 , e) 
+ c.c. 



The above expression can be rewritten using the magne- 
tization M defined below Eq. [571 



■con3 

Je i 



eJTfT 
' 2T< 



^7 T f da 

^ieij-^—cM z lbxi / —g<(p,e;A,imp) 

Jo P^°J 27T 

^- / ^dr 2 ...dr 4 n F (e)e ZQ/3 (r 2 + r 3 ) ce i; 3 (r 4 , r 4 ) 

o J 27T 



[.90* (ri , r 2 , e)crf 9 (r 2 , r 3 , e)gg (r 3 , r 4 , e) 
-9o(*u r 2, e)ed(r 2 ,r 3) e)#f (r 3 ,r 4) e) 



(B3) 



2T 7 (2^) 2 
x ^ (r 4 , rjgff (r 4 , r 2 , e)^f (r 3 , r 2 , we) gjf (r 3 , r 4 , e) 
x [L /? '(r 2 ,r 3 ,t i ;)n_ F (e - u>) - L A (r 2 ,r 3 ,ui)n F (e)] 
+ c.c. (B4) 

Here, we dropped terms with three retarded (advanced) 
quasi-particles Green functions and we used the identi- 
ties for the products of distribution functions presented 
in Eq.[50l We show now that 'f e em contributes to the cur- 
rent only at finite temperature, 'f e em ► 0. We demon- 
strate that jg em is closely related to the contribution of 
the Aslamazov-Larkin diagram to the magnetic suscep- 
tibility^ To calculate jg em , we allow the magnetic field 
to depend on the coordinate; in the end of the derivation 
we shall take the limit of a uniform magnetic field. 

Let us first concentrate on the product (r + 
r')L R ' A (r, r',u). Since the magnetic field varies in space, 
the dependence of the propagator L on the center of mass 
coordinate is through the magnetic field. The result of 
acting with the operator R a on L(R; p, ui) is equal to 
the derivative with respect to the vector potential Ap (r) 
multiplied by r a Ap(r): 



6 f dc f 

(r + r') a L R (r, r', uj) = — 4i- / — dr 1 dr 2 dr 3 ta,uh 

cj 2?r 



2T 



x L R (r, t u w)^(r 2 ,r 3 )r^(r 2 ) ff R (r 1 , r 2 , e')g R (r 3 , r 4 , e') 
xg A (r 1 ,r i ,Lu-e')L R (r i ,r',Lu). (B5) 



Now, we may return to the limit of a constant magnetic 
field, H(r) = Hz, and set the vector potential to be 
A(r) = H(r) x r/2. Inserting the expression given in 
Eq. IB5I into j^™ presented in Eq. IB41 we obtain: 



■rem _ ■ c v J ± 
Je i ~ lt V rp 



dede'duj 

— g— dr 2 ...dr 8 np(Lu) tanh 



(2tt) 3 



2T 



Szapv's (r s , ri)u^(r 4 , r 5 )r"e: 21//J #(r 4 )r 4 p (r 4 , r 2 , e) 



x SoVs, r 2 , u - e)g R {r 3 ,r i ,e) {L R (r 2 ,r 3 ,uj)g R (r 3 ,r 4 ,e')g R {r 5 ,r 6 ,e')g A (r 3 ,r 6 , u> - e r )L R (r 6 , r 7 , uj)n F (e - uj) 
+L A (r 3 ,r 2 ,uj)g A (r 3 , r 4 , e')g A (r 5 , r 6) e')g R (r 3 , r 6 , cj - e')L A (r 6 , r 7 , w)?i F (e)} + c.c. 

I 



(B6) 



The diagrammatic interpretation of the above expres- ducting fluctuations, the electrons are considered to be 
sion is presented in Fig. 1131( a) . Unlike the supercon- three dimensional. Then, for an isotropic system we 



22 



may rewrite the product of the antisymmetric tensors 
as: e za ps Z vn = (5 v , a 8n,0 - Sv,p5^ a )/3. Now we average 
over the disorder (see Fig. [I3jb)) and represent all the 
propagators (Green functions, propagators of the super- 



conducting fluctuations and Cooperons) as a product of 
the phases and the gauge invariant terms. All the phases 
are collected into the function e J *. After performing the 
Fourier transform, the equation for j^™ becomes: 



• ^V/T f dede'dm dkdk'dqdQ ( uj - e' \ i$ d 2 H(Q) _ fl 
T* 4 — — \0v aOu. b — o v rOu, a) I —, — ^ ; — 7T3 — np(w)tanh e — — <7n (k, e) 

x v$(k, k + Q)g R (k + Q, e)^(q - k, u - e) {c R (q + Q,e,w- e)L H (q + Q, wJC^q + Q, e', w - e') 

x.g i? (k' + Q, e'K(k' + Q, k')5 fl (k', e')ff A (q - k', w - e')C i? (q, e', w - e')i fl ; (q, w)C fl (q, e,w - e)n F (e - w) 
+C**(q + Q,e,w- e)L A (q + Q, ^)C* A (q + Q, e', - e')g A (k' + Q, e'K (k' + Q, k')g A {k' 7 e')g R (q e') 
xC A (q, e', w - e')L A (q, w)<7 fl (q, e, w - e)njr(e)} + c.c. (B7) 



Here, H(Q) = (27r) d iJ5(Q)z and v (k',k) = 
f^-i-xl + r-fi-x mr) is the Fourier 

\ 2m c ak 2m c ok.' J 

transform of the velocity; the arrows above the deriva- 
tives indicate on which of the Green functions the 
derivative acts. 

Let us transfer the derivatives with respect to the mo- 
mentum Q from the magnetic field to the rest of the 
expression using integration by parts. Since $ which col- 
lects all the phases contains only derivatives with respect 
to the momenta, it will not be differentiated in the course 



of this operation. In the diffusive limit, the main contri- 
bution to the current is obtained when the derivatives 
with respect to Q act on the propagators of the collec- 
tive modes, either L or C (here we rely on the argu- 
ments given below Eq. |44|) . Moreover, it follows from 
the tensor structure of Eq. IB7I that only terms in which 
d 2 / dQ a dQv = 26 a . u d/dQ 2 survive. [One should keep 
in mind that the gauge invariant propagators of the su- 
perconducting fluctuations and Cooperons depend on the 
square of the momentum.] Then, jg 6 ™ can be written as: 



■rem _ ., ■ 

Je i ~ ^'t'ij 



e 2 V,T f dede'duo dkdk'dq 



3T c 



-H 



(2tt) 3 (27r) 3d 
d 



np(uj) tanh 



2T 



e**5o«(k, e K(k,k) 



x g R (k,e)g A (q- k,w - e) \ ^ [C R (q, e, w - e)L R (q, w)C R {q, e', u - e')\ ~g R (k', e')v%(k', k').g fl (k', e') 



Xg A (q - k', w - e')C R (q, e',U- e)L R (q, uj)C R {q 7 e, CJ - e)n F (e - u 



_d_ 

dq 2 



C R (q, e,U - e)L A (q, uj)C A (q, e> - e>) 



xg A (k' 7 e'K(k', k')ff A (k', e')g R (q -k',co- e')C A (q, e' , w - e')L A {q, uj)C R (q, e,u- e)n F (e)} 



+ c.c. 



(B8) 



The next step is to integrate over the electronic de- 
grees of freedom and to transform to the basis of the 
Landau levels. This can be performed following the ex- 
planation presented in Sec. IIVI The only difference is in 
the matrix elements for the Landau levels. While in the 
calculation presented in the main text the matrix element 



is (N, 0\V x V y \M, 0), here we have (N, 0\V 2 + V 2 \M, 0) = 

AD 2 [(N + l)<W-i + (M + l)cW-i] /4- Finall y> af- 
ter replacing the derivative with respect to q 2 with a 
derivative with respect to the index of the Landau levels, 
the expression for j\ em acquires the form: 
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3lT = -4™^%/^ J £(JV + 1) tanh n» ± [n F {e - u )C$(e, . - e)L%( U ) 

xC§(e\u;-e')C% +1 (e',u-e')LX +1 (w)C% +1 (e,cu-e)+n^^ 

x Li +1 (w)C# +1 (e, u-e)]+ c.c. (B9) 



To recognize that the above expression goes to zero in 
the limit T — > 0, we have to integrate over the fermionic 
frequencies e and e': 



•rem 
Je % 



24tt 2 



2 eVj-T 



/ dw^(iV + l)rip(w) (BIO) 



iuo Q c (N + 3/2) 



1 

2 4ttT 

2 



icj fi c (AT + l/2) 



4^T 



4ttT 



4ttT 



+ c.c. 



now. we can exploit the fact that in the expression dif- 
ferentiated with respect to N the argument N stands to- 
gether with the frequency, — iuo + Q C N. Replacing the 
derivative d/dN by a derivative with respect to the fre- 
quency and integrating by parts, one immediately gets 
that the above integral includes the factor dnp(uj)/dui 
and, hence, vanishes at zero temperature. 

To summaries, we showed that in accordance with the 
third law of thermodynamics, the transverse thermoelec- 
tric transport coefficients a xy and a xy go to zero in the 
limit T — ► 0. [As to the longitudinal coefficients, both 
jcon3 an( j ynag j nQ ^ a pp ear w hile the remaining contri- 
butions vanish independently at T — ► 0.] It follows from 
this result that at finite T, one can obtain the tempera- 
ture dependence of the thermoelectric current by substi- 
tuting np(u>) with np{uj) + 0(— uj) in all the expressions 
determining j™" andj" 1 " 9 . 



APPENDIX C: THE QUASI-PARTICLES 
CONTRIBUTION TO THE THERMOELECTRIC 
TRANSPORT COEFFICIENTS 

In this Appendix we discuss the role of the particle-hole 
asymmetry and the constant density of states approx- 
imation in determining the quasi-particles contribution 
to the thermoelectric transport coefficients. In metallic 
conductors the quasi-particle excitations yield a negligi- 
ble contribution to the Ncrnst effect and to its counter- 
part, the Ettingshauscn effect. Let us consider a system 
with an electron and hole conducting bands that have a 
particle- hole symmetry, i.e., the bands are identical and 
the two species of particles differ only in their charge. We 
shall describe the deviation of the distribution functions 
of the two species (<5/ e and Sfh) from their equilibrium 
value in the linear response to a temperature gradient. 
For that we write the classical Boltzamann equation in 



the relaxation time approximation: 

SfeA e K) <9/ (e k ) n „ T ev k xH5/o(£ k ) 
-v k Vi =F — ■ 



dT 



dk 



(CI) 



Here the equilibrium (Fermi-Dirac) distribution function 
is denoted by /o(efc)j and Vfc is the velocity of the parti- 
cles. 

The electric current is the sum of the electric currents 
due to the electrons and the holes: 



itotal 



f dk f dk 

~ 2e J (2^p Vk<5/e(ek)+2e y (2^p Vk5//i(ek) 



(C2) 

Notice that the factor of 2 results form the sum over 
the two spin directions. For simplicity we only examine 
the limit of vanishingly small magnetic field. In order 
to determine whether a current vanishes in the particle- 
hole symmetric system we just need to count the powers 
of the electric charge; an odd power means cancellation 
of the two contributions to the current. 

We start from the longitudinal electric current induced 
by the temperature gradient. In the limit H — > 0, the 
longitudinal current is independent of the magnetic field: 



2c 



dk d/ (e k ) 



To 



[ekD e - e k D h ] = 0, 

(C3) 

where D e = D) x = D = v\r/d with d the dimension of the 
system. Since the expression includes only one power of 
the charge, there is no longitudinal electric current unless 
particle-hole asymmetry is introduced. 

The transverse current is obtained when the Lorentz 
force in the Boltzmann equation acts on the distribution 
function. Therefore, the expression for the transverse 
current contains an additional power of the charge: 



2e 



(2ir) d de k 



dk S/ (e k ) V X T 

-e k D [lu c t - {-w c t)\ ——— + 0. 

^0 

(C4) 

The additional charge enters through the cyclotron fre- 
quency lu c = eH/m*c. The even power of the charge 
means that the particle-hole symmetry does not con- 
strain the Nernst effect. 

Now, we look at the contribution for the transverse 
electric current in a metal with only one conducting band. 
We use the approximate of a constant density of states 
which is standard for Fermi liquid systems. Under this 
approximation the expression for the transverse current 
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is: 

j e « = 2eu Q D(cj c r)^ f de k ^^e k . (C5) 
To J <9e k 

Since near the Fermi energy the integrand is an odd func- 
tion of the energy, the resulting current is zero. There- 
fore, under the approximation of a constant density of 
states at the Fermi energy this contribution vanishes^ 
One may conclude that in metallic systems with high 
Fermi energy the contribution of the quasi-particles to 
the Nernst signal includes a small factor related to the 
deviation from the constant density of states which is of 
the order T/ep. In semi-metals like Bi where the con- 
stant density of states approximation cannot be used, a 
large Nernst signal was measured^ 

Let us compare the magnitudes of the transverse 



Peltier coefficient generated by the quasi-particles and 
by the superconducting fluctuations in a film of thick- 
ness a. The first is of the order ~ {oj c T)evDaT / £p for 
uj c t <C 1 while the second one is of the order ~ ett c /T 
for Q c /T -c 1 and ~ eT/£l c for the opposite limit. Thus 
in the limit of vanishing small magnetic field the ratio 
between the contribution of the quasi-particles and the 
fluctuations is a™ l/a{.y ~ (kpa)T 2 r/ 'ep. At higher mag- 
netic fields (but still in the limit oj c t <C 1) this ratio 
becomes a%y/al y ~ (kpa)epT(uj c T) 2 <C 1. Under the 
condition of the experiment^ the ratio a%P/al l y -C 1 
up to T < 100T C and H < 100H C2 . The reason why the 
Nernst signal generated by the superconducting fluctua- 
tions dominates the one produced by the quasi-particles 
was explained in the end of Sec. IIVI 
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